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COMPUTER PROGRAM FOR DETERMINING ROTATIONAL LINE INTENSITY FACTORS 


FOR DIATOMIC MOLECULES 
Ellis E. Whiting 
Ames Research Center 


SUMMARY 


A Fortran IV computer program, that provides a new research tool for 
determining reliable rotational line intensity factors (also known as Honl- 
London factors ) , for most electric and magnetic dipole allowed diatomic tran- 
sitions, is described in detail. This "users' manual" includes instructions 
for preparing the input data, a program listing, detailed flow charts, and 
three sample cases. The program is applicable to spin-allowed dipole transi- 
tions with either or both states intermediate between Hund's case (a) and 
Hund's case (b) coupling and to spin-forbidden dipole transitions with either 
or both states intermediate between Hund's case (c) and Hund's case (b) cou- 
pling. It is not applicable to quadrupole transitions or to transitions 
involving an electronic state approximated by Hund's case (d) coupling. 


INTRODUCTION 


This paper describes a comprehensive computer program for the determina- 
tion of rotational line intensity factors (also known as Honl-London factors) 
of diatomic molecules. The program is based on the theory of the intensity of 
rotational lines in diatomic molecular spectra presented in references 1 and 2 

The first analytic formulas for the rotational line intensity factors of 
diatomic molecules were derived by H. Honl and F. London (ref. 3) in 1925 for 
the simple singlet-singlet electronic transitions. Since that time satisfac- 
tory general formulas for the intensity factors of singlet, doublet, and 
triplet transitions have been developed. (See ref. 1 for a brief history of 
the development of analytic intensity factors.) 

The availability of realistic intensity factors for all spin-allowed or 
spin- forbidden transitions is highly desirable. 1 However, the algebraic dif- 
ficulties associated with deriving general formulas for intensity factors, 
when the multiplicity is four or greater, make it unlikely that tractable 
formulas will be derived for them in the foreseeable future. Therefore, the 
provision of numerically determined intensity factors appears to be the only 
practical and acceptable alternative. 


^or spin-allowed transitions AS = 0, AA = 0, ±1, and For 

spin- forbidden transitions Afi = 0, ±1 and at least one of the spin-allowed 
conditions is violated. 



The computer program described is a new research tool for determining 
reliable intensity factors for most electric and magnetic dipole-allowed 
diatomic transitions by numerical methods. It is applicable to spin-allowed 
dipole transitions with either or both states intermediate between Hund's 
case (a) and Hund's case (b) coupling and to spin-forbidden dipole transitions 
with either or both states intermediate between Hund's case (c) and Hund's 
case (b) coupling. The program is not applicable to quadrupole transitions or 
to transitions involving an electronic state approximated by Hund's case (d) . 

This paper is intended to provide an adequate users' manual for the com- 
puter program. Therefore, it includes a brief summary of the theory, a dis- 
cussion of several important program operations, three sample cases, detailed 
flow charts, and a complete program listing liberally annotated. 

The author would like to acknowledge the invaluable discussions with 
Dr. J. T. Hougen during the course of this work and the cheerful assistance of 
Cheryl Whiting in proofreading the detailed flow charts. 


THEORY 


The theory of rotational line intensity factors is derived from first 
principles, and in great detail, in reference 1. That derivation will not be 
duplicated here. However, in order to understand many of the definitions and 
operations used in the computer program, the user must be familiar with cer- 
tain key concepts in the derivation. Therefore, the following abbreviated 
discussion of the theory is included to provide the minimal information needed 
to use the program with some degree of confidence. 

Many of the terms used in the theory are defined in figure 1, a typical 
energy level diagram of an electronic transition from a 3 n to a 3 Z electronic 
state. The electronic spin-splitting and A-doubling, shown in the central 
portion of the figure, are greatly exaggerated. The figure illustrates the 
well-known fact that there are (25 + 1) electronic substates for sigma elec- 
tronic states and 2(25 + 1) electronic substates for nonsigma electronic 
states. The number of substates for both cases can be written in a unified 
fashion as (2 - 6 Q (25 + 1), using the Kronecker delta symbol. 

The structure of each electronic substate is composed of a series of 
vibrational levels each containing many rotational 5-levels. Further, each 
rotational 5-level can be split by a magnetic field into (25 + 1) Zeeman 
states. 

The basic spectral elements are Zeeman components produced by transitions 
between Zeeman states. The sum of all Zeeman components between two 5- levels 
•is a rotational line. This definition of a rotational line is identical to 
that of an atomic line (ref. 4, p. 237) and provides a natural bridge between 
atomic and molecular theory. It is also essential for the description of the 
theory in terms of basis functions discussed later. According to this defi- 
nition of a rotational line, a A-doublet is composed of two rotational lines. 
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However, in the literature a rotational line is frequently defined to include 
both components of a A-doublet (ref. 5 and 6 for example) . This difference 
provides one of the important sources of confusion in the theory of diatomic 
spectra. 

A vibrational band is the composite of all rotational lines occurring 
between vibrational levels v' and v" in all electronic substates. Figure 1 
shows only that part of a vibrational band formed between two substates. 

A band system is the composite of all rotational lines between two electronic 
states. Thus, a band system in molecular spectra is equivalent to a very 
extensive atomic multiplet in atomic spectra. 

The effect of nuclear spin is not explicitly included in the analysis. 
However, in nearly all cases, the nuclear spin and the resultant angular 
momentum without nuclear spin commute. Thus, the nuclear hyperfine line com- 
ponents can be obtained from the rotational line intensities by the methods 
of angular momentum theory summarized in reference 4. 

The theory of rotational line intensities is closely dependent on the 
various angular momentum operators and the way they couple to form their 
resultant. The vector model coupling diagrams for Hund's cases (a), (b) , and 
(c) coupling and their appropriate selection rules are shown in figure 2. In 
this discussion, Hund's case (a) basis functions are used for spin-allowed 
transitions and Hund's case (c) basis functions are used for spin- forbidden 
transitions . 


The starting point for all discussions of spectral intensities is 
Einstein's phenomenological equation for spontaneous emission in either an 
atomic line or a molecular rotational line (ref. 7) ; that is 


E , i T n = N T ihv T ! ,h 
d d d d d d d 


where 

the emitted power/unit volume, W m -3 

Nj t the population density of the j' rotational level, m -3 

h\)j\j\\ the energy of each photon emitted, W s 

the spontaneous emission transition rate per particle, also called 
the Einstein A coefficient, s _1 


In reference 8 Dirac derived the Einstein A coefficient in terms of an 
expansion in the electric and magnetic moments within an atom or a molecule. 
The strongest term in the expansion, if it is nonzero, is that due to the 
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( 2 ) 


electric 2 dipole moment, P; that is. 


Vj" 



(2 J' + 1) 


(3) 


where 

^ , the wave function of the upper Zeeman state M' 

the wave function of the lower Zeeman state A/' 

? the electric dipole moment operator in the laboratory coordinate 

system, C m 

e Q the permittivity of a vacuum, 8. 854xlO -12 C 2 J“ 1 m -1 

S T , ,m the line strength, defined in reference 4 

d d 


With equation (3), equation (1), for the power emitted in a line by 
spontaneous emission, can be written 


E J'J" 


3 4 

16lT Vj" 


V 


S J'J" 


(2 J' + 1) 


(4) 


The advantage of using equation (4) instead of equation (1) is that the line 
strength Sj'j" is symmetrical in the upper and lower states. Thus, the line 
strength is the same in emission and absorption (see ref. 4, p. 98). 

In order to make further analytic progress in the theory of diatomic 
molecules, the Bom-Oppenheimer approximation is introduced (see ref. 9). 

This approximation assumes that the total wave function can be written as a 
product of an electronic wave function and a vibrational-rotational wave 
function. Further, Pauling and Wilson (ref. 10) have shown that the 
vibrational-rotational wave function can also be approximated by a product of 
vibrational and rotational wave functions. On the basis of these approxima- 
tions, it is shown in reference 1 that the line strength for isotropic radia- 
tion can be written 


= (5) 


2 The derivation given herein and the computer program are based on the 
electric dipole moment. However, the results can also be applied to magnetic 
dipole transitions (see p. 19). 
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where 


<7 , ,, the vibrational Franck-Condon factor, dimensionless 

R the electronic transition moment, C m 

e 

a the rotational matrix element including the summation over M' and M" , 

dimensionless 

ip jpty-r the electronic-rotational wave functions for the rotational levels 

J' and J" 

R e ‘7!%* V 

for perpendicular transitions, that is, those for which An = ±1. 



for parallel transitions, that is, those for which Aft = 0, where p , p , and 

p are the components of ? in the molecular coordinate system. x ^ 
z 

The matrix brackets applied to Sjiju in equation (5) indicate that in 
diatomic spectra there may be several rotational lines in the same vibrational 
band with the same values of J ’ and J" (see fig. 1). Further, the symbolism 
correctly implies that the strength of the lines can be found by regular 
matrix operations. 

In general, it is not possible to write simple expressions for the matrix 
elements of <R e a> between the wave functions ipy and \l>£. Therefore, the 
solution proceeds in three major steps. First the wave functions are expanded 
in terms of a complete set of simplified basis functions $ (see ref. 2); next 
the <i? e a> matrix is expressed in terms of the basis functions; and finally 
the <i? e a> basis matrix is transformed to the wave functions. 

The transformation from the <R e a> basis matrix to the <l? e a> matrix 
between wave functions is given by the unitary transformations that symmetrize 
the basis functions ^sym and that diagonalize the Hamiltonian matrix between 
symmetrized basis functions Tfc a g. Thus, 

<*uK a K> s (6) 

where 

$ = |A£T>|fiJ> 

is the product of an electronic and a rotational basis function with Q = A + E 
and T = ^ S ym T di a g’ 
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The electronic transition moments R e are usually determined by comparison 
of theory with experimental data; they are, therefore, carried in the present 
analysis as unknown parameters. To calculate realistic line intensities the 
user must know these parameters. 

The rotational matrix elements, a, are related to the direction cosine 
transformation, which expresses the components of the dipole moment operator 
in the laboratory coordinate system in terms of its components in the molec- 
ular coordinate system. The matrix elements of the direction cosines are 
expressed in terms of the rotational quantum numbers in reference 11, and 
these are used in reference 1 to give the expressions for a (valid for 
isotropic radiation) listed in table 1. 

The general selection rules A J = 0, ±1 and A £2 = 0, ±1 are explicitly 
indicated in table 1. Separate matrices for the P-, Q - , and P-branches are 
formed, depending on the value of A J. The specific matrix elements for each 
case are given by the expressions in the appropriate column of table 1. 

An example <P e ct> basis matrix for a spin-forbidden 2 II ■*- 4 £ + transition 
is shown in table 2. The rows of the matrix are designated by the final state 
basis functions and the columns by the initial state basis functions. The 
number of independent matrix elements can be reduced because symmetric matrix 
elements are related. Symmetric matrix elements are those that interchange 
when the signs of both A and 2 are changed in both the upper and lower elec- 
tronic states (Rq and R e a 1Q in table 2, for example). 

110 ^ 

The relationship between the symmetric matrix elements is determined by 
the phase factors that occur naturally in the quantum mechanical description 
of the problem. For phase conventions consistent with those chosen by Condon 
and Shortley in reference 4 and by Hougen in reference 2 (eq. 2-31), the 
following "symmetry rules" are derived in reference 1: 

1. All the transition moments are real or all are pure imaginary and, 
for all practical purposes, they can always be chosen to be real. 

2. All of the symmetric matrix elements in any given matrix are either 
equal to or the negatives of each other. 

3. The symmetry in the P- and P-branch matrices is always opposite to 
the symmetry in the ^-branch matrix. 

4. The symmetric matrix elements in the ^-branch matrix can always be 
chosen to be equal, except for the following two cases when they must be 
chosen to be negatives: 

+ + , 

E •«-»-£ transitions with AS = 0, 2, etc. 

E*-*-^* transitions with A S = 1, 3, etc. 

As an illustration of these rules, table 2 is recast in table 3, as it would 
appear for a P-branch matrix. 


6 



In the case of spin-allowed transitions, the additional selection rule 
A£ = 0 applies to R e between electronic basis functions. This additional 
restriction reduces the number of nonzero matrix elements in the <f? g ct> basis 
matrix for spin-allowed transitions to those along the diagonal, where 
An = AA, as n = A + E. 

Furthermore, for spin-allowed transitions it can be shown (see ref. 2) 
that the absolute values of all of the nonzero electronic transition moments 
Rq are equal. Thus, for spin-allowed transitions, R e can always be factored 
out of the <R 0 a> basis matrix, and the line strength (eq. (5)) can be 
written as 

and the matrix elements of <<^ j'j" > > which are the rotational line inten- 
sity factors or Honl- London factors, only involve the transformed rotational 


matrix elements a; that is. 





(9) 

= 


(10) 


For spin- forbidden transitions it is not possible, in general, to sepa- 
rate the electronic and rotational parts of the problem. For these cases 
reference 1 shows that the transformed < R e a > matrix elements can be 
written 


= <^“1 + \ a 2 + • * C > (11) 

where cT. is the transformed rotational matrix element associated with R 

V 

Equation (11) illustrates the property that for spin-forbidden transi- 
tions the matrix elements of < R e a > cannot, in general, be separated into 
rotational and nonrotational factors. However, if one of the transition 
moments (denoted by R e ) is factored out of the <R e a> matrix, the expres- 

sion for the rotational line strength is similar to that for spin-allowed 
transitions; that is, 

< 5 JV"> = Vv’-^K + + ™k + * * * >! 2 (12) 

U 
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where 


P 




R 

e k 

T ■ g- * etc - 

e i 


The factors a, t, etc., are called intensity parameters. The squares of the 
matrix elements in equation (12) are called the rotational line intensity fac- 
tors and are defined as 


= I < PO£ + dS. * ™ k + * * ‘ > | (13) 

If equations (12) and (13) are combined, the rotational line strength for 
spin-forbidden transitions can be written in the same form as spin-allowed 
transitions (eq. (8)); that is, 

= q v^v'' R e l < ^ 

The sum of the intensity factors for a given value of J' or J" is: 


1. Spin-allowed transition 

tf ) = (2 - 6 0)(a , +am) )(2S + 1)(2J + 1) (15) 

where <S , = 1 if a = fc and = 0 if a t b. 
a 3 b 

2. For each transition moment in spin-forbidden transitions 

= C(2J + 1) (16) 

where C - 1 if the transition moment occurs only once in the < R e a > basis 
matrix and C - 2 if the transition moment occurs twice in the < R e a > basis 
matrix. The only transition moments that occur only once in any < R e a > 
basis matrix are those for the O' = 0 to = 0 transition of £<-+£ transitions 
of odd multiplicity. This transition moment is nonzero for transitions 

with AS = 0, 2, etc. and for E ± -*-»-E + transitions with A S = 1, 3, etc. The 
theory is explained in more detail in references 1 and 2. 


AUXILIARY PROGRAMMING CONCEPTS AND DETAILS 


The major concepts involved in the program are illustrated by the concep- 
tual flow chart shown in table 4. In brief, the computer solution proceeds in 
the following logical steps. 
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1. Set up the upper and lower Hamiltonian matrices for each value of J' 
and J" . 

2. Symmetrize and diagonalize the upper and lower state Hamiltonians. 

3. Set up the relevant rotational matrix for each pair of rotational 
levels J ' and J" . 

4. Transform the rotational matrices with the same transformations that 
transformed the upper and lower state Hamiltonians. 

A complete listing of the computer program is given in appendix A. The 
listing is liberally annotated, and if it is read in conjunction with the 
flow charts in appendix B, it should be nearly self-explanatory. In any com- 
plicated computer program, however, there are always a few programming details 
for which the logic is not immediately obvious. The following topics are 
included to give the reader some insight into the more obscure details. 


Initial State to Final State Notation 

The computer program was organized at its inception to describe the tran- 
sitions from the initial state (columns of the matrices) to the final state 

(rows of the matrices) . Unfortunately, this choice complicates the logic 
necessary to print correct titles for the calculated results. However, the 
program works, and the substantial changes necessary to switch to the more 
standard notation, which describes the rows of the matrices with the upper 

state and the columns of the matrices with the lower state, does not seem jus- 

tifiable at this time. Further, the rows and columns of all basis matrices 
are ordered from top to bottom and from left to right in terms of the basis 
functions |AS , Z>|jU> as follows: 

|+A S +S>\A+S J > , |+A S +S-1 >|A+S-1 J > , • • • , |+A S -S )>|A-S J >, 
|-A 5 +5>|-A+5 J>, | -A S +S-1 > | -A+S-l J>, • • • , | -A S -S > | -A-S J > 

Absorption and/or Emission 

The intensity factors can be calculated for either absorption (in terms 
of J") or emission (in terms of J '). The accepted standard notation for rota- 
tional lines is always J" . However, since in quantitative calculations 
involving line emission it is usually more convenient to denote intensity 
factors by J' , this flexibility was provided in the program logic. 


Hamiltonian Matrix 

Hamiltonian or energy operator matrices are set up for each value of J in 
both the initial and final electronic states. That is, the energy levels of 
all rotational levels with the same value of J , in a given vibrational level 
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and electronic state, are collected into a single matrix. In general, when 
J > A + S, the Hamiltonian matrix contains (2 - 6 0 ^) (25 + 1) rows and col- 
umns. However, if A ^ 0 the two submatrices for +A and -A are mirror images. 
An example case is shown in table 5 for a 2 n electronic state. As there are 
no off diagonal terms between +A and -A, these two submatrices do not inter- 
act. It is therefore only necessary to operate on one of these submatrices, 
and, when needed, the full operator matrix can easily be constructed. 

The two operations performed on the Hamiltonian matrix between basis 
functions are those that symmetrize the basis functions and that diagonalize 
the Hamiltonian matrix between symmetrized basis functions. The Hamiltonian 
matrices between basis functions and between symmetrized basis functions are 
identical for nonsigma electronic states. The unitary transformations that 
transform the Hamiltonian matrix to symmetrized basis functions for 3 E and 4 E 
electronic states are shown in figure 3. The generalization to any multi- 
plicity is straightforward. 

The diagonalization of the Hamiltonian is performed by the EIGEN subrou- 
tine, which finds T diag for symmetrical, real matrices. This subroutine is a 
slightly modified form of the EIGEN subroutine described in the IBM System 360 
Scientific Subroutine Package, document H20-0205. 

In this program we are concerned only with determining the rotational 
intensity factors and do not solve explicitly for the energies of the rota- 
tional levels. Therefore, to a very good order of approximation, it is only 
necessary to include in the Hamiltonian the major energy interaction terms. 

For nonsigma states only the first-order spin-orbit interaction term, AM, is 
included. For sigma states, the first-order spin-orbit interaction term is 
zero; therefore, both the second-order spin-orbit term and the spin-spin 
interaction term are included. Both terms produce a similar effect and are 
lumped together as A E, the energy separation between spin states extrapolated 
to J = 0 (see appendix C) . If A (or A E) is negative the state is an inverted 
state. 

The unitary transformation matrix that diagonalizes the Hamiltonian is 
not affected by a constant value along the diagonal of the Hamiltonian matrix 
or by a constant times every matrix element. Therefore, as we do not need the 
rotational energies themselves, constant or ^-dependent only terms along the 
diagonal are removed from the Hamiltonian and all matrix elements are divided 
by Bfi 2 . 

General expressions for the Hamiltonian matrix elements are given in 
reference 2. However, on the basis of the above discussion, only the follow- 
ing terms axe included in the computer program. 

1. Nonsigma states 

H(K,K ) = -fi 2 - E 2 + YM 

H(K,K+ 1) = -[(el - + Q + 1 )(S - E)(S + E + l)] 1 / 2 

H(K+l,K) = E{K,K* 1) 


10 


(17) 

(18) 
(19) 



2. Sigma states 


H{K,K) = -fi 2 - E 2 + A E/B (20) 

#(F,F+1) and H(K+l,K) same as above 


where Y = A/B, and K specifies the row and column indices with (1,1) as the 
upper left matrix element. 

After the Hamiltonian matrix is diagonalized, the largest energy level 
is in the upper left matrix element and the smallest is in the lower right 
matrix element. This organization specifies the final form of T^^g and, 

' hence, determines the order in which the branch lines occur in the intensity 
factor matrices. In order to understand this point, it is necessary to know 
how the rotational levels and the branches are designated. 


Designation of Rotational Levels 

In the standard notational scheme, the rotational levels of diatomic 
molecules are designated by two parameters in addition to J : (a) the rota- 
tional quantum number N , exclusive of spin and (b) the spin substate 
^ 1 ,^ 2 , . . . , F 25 + 1 . The designation of rotational levels by N is most 
j appropriate for Hund's case (b) coupling, where If is a valid quantum number, 
j However, there is a one-to-one correspondence between the rotational levels in 
Hund's case (b) and any other coupling case, so that a value of N can always 
be assigned. 


The designations of N and F. for the spin substates, when J > A + S, are 
related as follows: % 


F 1 

=> 

S3 

II 

+ s 

or 

52 

II 

l 

F 2 

• 

=> 

• 

II 

52 

+ 5-1 

or 

I 

II 

83 • 


( 21 ) 


25+1 


=*> J = N - S 


or N = J + S 


From these equations it is clear that for the group of F ^ levels with the 
same value of J, F± corresponds to the lowest value of N and F 2 g + i corresponds 
to the largest value of N. Because N = ri + R, the difference between these 
N values is due to the difference in rotational energey, R. Thus, for the 
group of F^ levels with the same value of J, the rotational energy increases 
from fj, the lowest energy level, to ^S+l’ the highest energy level. 


In nonsigma electronic states (A ^ 0), the assignment of D, where 
fi = A + E, to the F £ levels depends on whether the electronic state is regular 
or inverted. In regular states the smallest value of fi is associated with the 
lowest energy level (i.e., Fj) . In inverted states the opposite associa- 
tion is made. Therefore, the assignment of N , F . , and J to the rotational 

ts 
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levels of nonsigma electronic states can be made in the following empirical 
fashion. 


1. For J > A + S and for rotational levels with the same value of J, F 1 
is assigned to the lowest energy level, F 2 to the next higher energy level, 
etc. 

2. The value of ft is assigned to each F. substate based on whether the 
electronic state is regular or inverted. 

3. The rotational quantum number J is assigned sequentially from the 
lowest rotational level in each F. substate, where the minimum value of 

J is |n| . 


4. N is assigned as specified in equations (21) with the restriction 
that N > A. 

The steps outlined above are applied to the rotational levels of a 4 n 
electronic state in figure 4. The value of N is shown to the left of each 
row, and the value of J is shown on the line representing the rotational 
level. The separations of the energy levels are not drawn to any physical 
scale, but they do indicate that for a given value of J, the energy increases 
with increasing N. 

A study of figure 4 shows that, for either regular or inverted 4 n elec- 
tronic states, there are only two rotational levels with J = 1/2 and only 
three rotational levels with J = 3/2. The full multiplicity is therefore not 
developed until J > A + S. Also, note that the N designation of the lower 
rotational levels in regular states is not given. The empirical scheme breaks 
down for these levels. To find the appropriate value of N, when J < A + S, 
one must operate on the wave function with N 2 . The substate designation, F., 
is established for J > A + S and is extended to low J levels as described 
above . 

Much of the discussion above for nonsigma states also holds for sigma 
states (A = 0) , but the designation of the rotational levels with N, F and 
J is the same in both regular and inverted states. The assignment for a 
3 E electronic state is illustrated in figure 5. The concept of ft used for 
nonsigma states is not valid for sigma states and is not shown. 

A study of figure 5 shows that there is only one rotational level with 
J = 0 and that it fits naturally into an assignment of F 3 . This assignment of 
J = 0 agrees with that shown by Herzberg (ref. 12, p. 223) and also is compat- 
ible with Hougen's assignment of F 3 and FV to the J = 1/2 levels of 4 E states 
(ref. 13). However, Tatum and Watson (ref. 14) chose to assign the J = 0 
level of 3 Z states to F j for regular states and to F 3 for inverted states. 
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Designation of Branches 


The standard scheme used to designate the branches in each vibrational 
band, is based on using letters to indicate the changes in J and N occurring 
during the transition and on including the F ^ assignment of the upper and 
lower substates. The assignment of letters to indicate the values of A J and 
AW is summarized in figure 6. The selection rules for dipole radiation limits 
W to 0 or ±1. The branch designation scheme is illustrated in figure 7. If 
AW = AJ, the upper letter is not included and if F\ = F \ !, only one subscript 
number is included (i.e., R R 2 2 = R 2 in fig. 7). 

The energy change during a transition is usually more closely associated 
with AW than with A J. Therefore, the form of the branch is also primarily 
controlled by AW. Hence, the R Q 32 branch is called the 7?-form ^-branch; that 
is, even though it is a Q-branch (A J =0), it usually has the form or appear- 
ance of an R - branch as AW = +1. 

The proper designation of the branches can be determined by forming 
matrices of the branches for fixed values of J ' , J" , and bJ . These matrices 
are illustrated in figure 8 for a spin-allowed 4 n«->- £f £ transition and in fig- 
ure 9 for a spin-forbidden 2 IH-»- 4 r[ transition. The rows of the matrices are 
labeled by the assignments for a given value of J' and similarly the col- 
umns by F V for a given value of J" . 

The designation scheme illustrated in figure 9, for the branches of spin- 
forbidden transitions, is not universally applied. For example, in refer- 
ence 15 Kovacs designates the F ^ levels in the 2 n electronic state in 2 IR-* 4 n 
transitions as F 2 and F 3 sublevels rather than F 1 and F 2> as shown in fig- 
ure 9. The designations F 1 and F 2 for the 2 n state, however, are consistent 
with the recommendation made by Mulliken (ref. 16), and it seems desirable to 
retain this designation for all types of transitions. 

The physical characteristics of rotational lines, such as their wave- 
lengths and their intensity factors, are not, of course, dependent on the 
notation used. Therefore, the designation of branches can be altered to suit 
personal preference by making appropriate substitutions in the branch symbols. 
This option also applies to the designation of the J = 0 rotational level in 
3 E states mentioned previously. However, standard designation schemes are 
very desirable. 


The Rotational Matrices 

The rotational matrices, in terms of the basis functions, are constructed 
from the matrix elements given in table 1. In the case of spin-allowed tran- 
sitions, the electronic transition moments are factored out of the <W e a> 
basis matrix (see eq. (7)), and all the matrix elements are transformed simul- 
taneously. Further, for spin-allowed transitions the rotational matrix 
elements are determined individually, rather than by the symmetry rules given 
on page 6. Also, for spin-allowed transitions the rotational line intensity 
factors are unchanged when the symmetrizing transformation is neglected and 
when only the +A submatrix is explicitly considered if A > 0. 
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In the case of spin-forbidden transitions, each rotational matrix element 
and its symmetrical counterpart are associated with a specific independent 
transition moment. Thus, each pair of symmetric matrix elements must be 
transformed separately by both the symmetrizing and diagonalizing transforma- 
tions (see eq. (6)) . Because only two matrix elements are involved in each 
transformation, the symmetrizing transformation in the program is performed 
algebraically rather than by the complete unitary transformation. 

Each matrix element in the transformed rotational matrices is associated 
with a particular, branch, based only upon its location within the matrix (see 
figs. 8 and 9). However, for transitions involving nonsigma electronic 
states, two submatrices occur and the usual branch designation scheme, 
described previously, does not at first appear adequate. However, if we print 
only one line intensity factor for each A-doublet in X++Y transitions (Z and Y 
represent any nonsigma electronic states), and if we combine the two sub- 
matrices in Z++X transitions prior to printing the results, the designation 
scheme described in the previous section is adequate. 

For spin-allowed transitions, the above simplification was introduced 
indirectly by including only the <ot> basis matrix elements from the +A 
submatrix and by neglecting the symmetrizing transformation. Thus, for spin- 
allowed transitions, all the required matrix elements occur in the upper left 
submatrix. For 2>+il transitions, however, A-doubling does not occur and, as 
only one submatrix of the electronic state is included, the matrix elements 
must be multiplied by Jl. 

For spin- forbidden transitions the situation is slightly more complicated. 
For example, transformed rotational matrices for 3 Z 3 Z transitions are illus- 
trated in figure 10. The X in these matrices represents the only possible 
nonzero matrix element and the F £ designates the rows and columns assigned as 
discussed on page 11. Clearly, in either of these matrices if the elements 
of the lower submatrix are added to the elements of the upper submatrix with 
the same values of F. ^ and F £, we will always add a nonzero to a zero value or 
vice versa. Furthermore, the resultant upper submatrix will contain matrix 
elements for all the branches. The discussion of 3 X -> 3 I transitions is simi- 
lar except that left and right submatrices replace upper and lower submatrices. 
The generalization to any multiplicity is straightforward. 

For 3 X++ 3 Y transitions (neither electronic state is a sigma state), the 
transformed rotational matrices are illustrated in figure 11. The two matrix 
elements for each A-doublet occur in symmetrical locations with respect to the 
center of the matrix. As we print only one component of a A-doublet, we can 
always place the required matrix elements in the upper left submatrix by add- 
ing the elements in the lower left submatrix to the elements in the upper left 
submatrix with the same values of F\ and F ". . 

t- ^ 

Therefore, for all cases the transformed rotational matrix elements are 
organized into the upper left submatrix for printing. For spin-allowed tran- 
sitions these matrix elements are squared before being printed to form the 
rotational line intensity factors, or Honl-London factors. For spin-forbidden 
transitions these matrix elements are the cT. values shown in equation (13) and 
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they must be printed without being squared. Before the output is printed, 
however, when J < A + S, the matrix elements frequently need to be shifted to 
ensure proper labeling. 


Shifting of Rotational Matrix Elements When J < A + S 

At the conclusion of the transformation operations discussed in the pre- 
vious section, the transformed rotational matrix elements are located in the 
upper leftmost portion of the matrix. The output section of the program 
prints branch headings corresponding to the location of each matrix element in 
the fully developed matrix. However, when J < A + S, the matrix is not fully 
developed and the matrix elements may not be in the proper positions to cor- 
respond to the headings that are printed. In these cases the matrix elements 
are shifted before being stored in the SAVE array. The SAVE array is used for 
temporary storage, before printing, and is discussed later. 

The reason for shifting the matrix elements and the logic employed for 
shifting is indicated in figure 12. This figure shows the J = 1 rotational 
levels in the +A submatrix of the Hamiltonian matrix for a 7 n electronic 
state. 

The F- level designations of the rows and columns of a fully developed 7 n 
matrix (i.e., J > A + S = 4) are shown in figure 12(a). As discussed on 
page 11, the highest energy level (F 7 ) is in the upper left matrix element. 

The program is written so that the three energy levels for J = 1 (£3, E 2 , and 
Fj) naturally occur in the upper leftmost portion of the matrix, which corres- 
ponds to the F 7 , F g , and F 5 levels. However, if we determine the proper 
designations of the J = 1 rotational levels, we see that in a regular elec- 
tronic state these levels should be designated F 4 , F 3 , and F 2 , and in an 
inverted electronic state, they should be designated F g , F 5 , and F 4 . Thus, to 
correspond to the designations of the fully developed matrix, the matrix ele- 
ments must be shifted three spaces for a regular electronic state and one 
space for an inverted electronic state. These shifts are indicated by the 
heavy lines in figure 12(a). 

The shifting of the rotational matrix elements when J < A + S is as 
follows: The columns are shifted by the shift necessary in the initial state 

Hamiltonian matrix and the rows are shifted by the shift necessary in the 
final state Hamiltonian matrix. Generalization of the logic described above 
to other electronic states is tedious but straightforward. 


Designation of Transition Moments in Spin-Forbidden Transitions 

The discussion leading to equation (11) shows that several independent 
transition moments may be present in spin-forbidden transitions. Each of 
these transition moments is explicitly identified in the computer program, but 
the designation scheme, contrary to the matrix operations, must be specified 
in terms of upper and lower states. The designation scheme used for this pur- 
pose is illustrated in table 6. 
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Table 6 shows the <i? e a> basis matrix for the ^-branches of a 2 n-*-»- 4 n 
transition. The symmetry rules given on page 6 have been used to equate the 
symmetrical matrix elements. As noted above, the rows of the matrix refer to 
the upper ( 2 n) electronic state and the columns to the lower (^n) electronic 
state. The opposite choice could have been made, but the choice made corres- 
ponds to conventional matrix nomenclature. If, instead, the rows of the 
matrix are designated by the final electronic state and the columns by the 
initial electronic state, as in the program matrix operations, the same tran- 
sition moment may have one designation in absorption and a different one in 
emission. This situation is, of course, not acceptable. Thus, the logic in 
the program at this point is rather complicated. 

The matrix elements of the <P g a> matrix in the basis functions and, 
hence, the transition moments R e are designated with two single digit numbers. 
The first number specifies the upper electronic substate and the second number 
specifies the lower electronic substate on which the transition moment oper- 
ates (i.e., i? e (upper, lower)). 

In terms of upper and lower electronic substates, the columns of the 
< R e a > matrix are numbered from 1, beginning with the rightmost column if 
A" = 0 and with the rightmost column in the left half of the matrix if A" # 0, 
as in table 6. Note the circled numbers at the top of the columns in the left 
half of the <i? e a> matrix in table 6. Similarly, the rows of the <P e a> 
matrix are numbered from 1 upward, beginning at the bottom row if A* = 0. If 
A' ^ 0, the rows are numbered from 1 upward in the top half of the matrix, 
beginning at the dividing point, and from -1 downward in the lower half of the 
matrix. Note the circled numbers at the left of the rows in table 6. 

The rotational line intensity factors of the two lines forming a 
A-doublet in spin-forbidden transitions are not equal if the line strengths 
contain at least one <P e a> matrix element between basis functions with oppo- 
site signs on A' and A" (i.e., <±A' |i? e ot | *A"» and at least one matrix ele- 
ment with the same signs on A' and A". The <±A' |i? e a|+A"> matrix elements 
of nonsigma to nonsigma spin-forbidden transitions are indicated in the desig- 
nation scheme by a negative number in the upper (i.e., first) index location 
of the transition moment. There are three such matrix elements in table 6: 
V I’ 2 ), and i? e (-2, 1) . 

The P-branch matrix of a i+ E H V->- 1+ E - transition, shown in table 7, has five 
independent transition moments. 


SAVE, ITRANI, ITRANF, NTRANR, NTRANQ, and NTRANP Matrices 

The rotational line intensity factors (Honl-London factors) for spin- 
allowed transitions and the transformed rotational matrix elements for spin- 
forbidden transitions are stored in the SAVE array until the calculation is 
complete. The SAVE array is three-dimensional and can be viewed conveniently, 
as shown in figure 13. The three dimensions correspond to (1) the number of 
independent transition moments in the transition (maximum of nine) , (2) the 
number of branches in a vibrational band (maximum of 150, there are 
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147 branches in spin-allowed septet-septet transitions), and (3) the number of 
rotational lines in each branch included in the calculation (maximum of 200) . 
In the dimension for storing branches, the i?-branches are stored in the first 
50 locations, the ^-branches in the second 50, and the P-branches in the last 
50. Any dimensions of the SAVE array can be changed to any desired value, 
limited only by the size of the computer memory. 

In the case of spin-forbidden transitions, the designations of the inde- 
pendent transition moments must also be stored. These correspond to each 
occupied row in the short dimension of the SAVE array for each value of J. 

The upper substate (or first) designation of the transition moments in the 
< PgOt > basis matrix is stored in the ITRANI array and the lower substate (or 
second) designation is stored in the ITRANF array. The number of independent 
transition moments for each value of J in the P-, Q - , and P-branches is stored 
in the NTRANR, NTRANQ, and NTRANP arrays . 


Input Data and Sample Cases 

This section is intended to give the user a general picture of the ease 
of operation, broad generality, and potential applications of the computer 
program. 


The data needed to initiate a calculation are: 

1. The resultant spin S of each state. 

2. The A value of each state. 

3. The ± symmetry for E states. 

4. The spin-orbit and/or spin-spin parameters (see appendix C) . 

5. The values of J . and J 

rmn max 

6. The type of calculation (i.e., emission or absorption). 

The format for the data input cards is illustrated in table 8. 

The computer output format is demonstrated by partial listings of the 
printed output for three sample cases in tables 9, 10, and 11. The sample 
case in table 9 is a spin-allowed 3 E + ■* 3 n transition; that in table 10 is a 
spin- forbidden 3 II 5 E“ transition; and that in table 11 is a spin- forbidden 
2 n 4 n transition. The information in the tables is discussed below. 

The heading at the top of table 9 indicates that it is an allowed 
3 E -* 3 n transition. The energy separation of the upper (E) state AE/B is 
specified as -10; therefore, it is an inverted state. The spin-orbit coupling 
(T = A/B ) of the lower (II) state is specified as 100; therefore, it is a 
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regular state. The paragraph in the heading refers to the lack of a univer- 
sally accepted convention for designating the low J levels of l states (i.e., 
el < 1 in this case) . As noted on page 12, for this case the J = 0 rotational 
level in the £ state is designated as an F 3 level. 

The values tabulated for each branch are the rotational line intensity 
factors (also called Honl-London factors) from Jrrtin = 0 to J max = 4 as speci- 
fied, in the input data. The column titled SUM is the sum of all Honl-London 
factors printed for a given value of J and for the F level designation 
repeated in the column headings. Therefore, SUM contains all Honl-London fac- 
tors in the R - , Q - , and P-branches from a given J level. Each value of SUM is 
printed three times, once in each P-, Q-, and P-subsection. In table 9, for 
example, the first (upper level) F designation is repeated in the column head- 
ings because it is an emission calculation. Therefore, for J = 2 and for the 
F^ rotational level, SUM is given by 

SUM = PI + QR12 + PP13 + Q1 + PQ12 + 0Q13 + PI + 0P12 + NP13 

This value of SUM is printed in the appropriate place in the SUM column for 
the P-, Q- , and P-branches. The Honl-London factors printed are for individ- 
ual rotational lines. Therefore, if A-doubling occurs in the spectrum, the 
value printed is for only one component. Thus, the value of SUM, which is 
only the sum of the Honl-London factors printed, is (2c T + 1) for £-*-»-£ and 
for nonsigma to nonsigma transitions, and 2(2J + 1) for £-<->n or IH-»-E 
transitions. 

Two spin- forbidden sample cases are included in tables 10 and 11 to 
illustrate results with and without A-doubling in the spectrum. The informa- 
tion printed in the headings for these cases is similar to that described for 
table 9 except that table 11 also contains a comment about A-doubling. The 
primary differences from the printout for the spin-allowed cases are in the 
tabulations. Normally, in spin-forbidden transitions there will be more than 
one independent transition moment, and the transformed rotational matrix ele- 
ments must be printed for each of them. 

The values tabulated for each branch are the transformed rotational 
matrix elements, cT^. The intensity factors are formed from these numbers as 
shown in equation (13), which is repeated here for convenience: 

= I < P“£ + + + ' * * > I 2 

where p, a, x, etc,, are the intensity parameters defined following equa- 
tion (12) . 

The PARTIAL SUM listed for the forbidden transitions is the sum of the 
squares of the cT^ values printed in the row to the right of the transition 
moment designation. The total sum for each transition moment is found by add- 
ing the partial sums from all branches containing that transition moment. 
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Lambda doubling occurs in the spectrum of the sample case given in 
table 11. However, the strengths of the two lines composing the A -doublets 
are not equal in all branches because three of the transition moments occur 
between basis functions with opposite signs on A' and A" (see page 16). These 
transition moments are designated byi? e (-2,l), R e (- 1,1), and R e (-1,2) in 
table 11. The intensity factor for the A -component not shown in the computer 
printout is obtained by changing the sign of these terms in equation (13) . 
In general, when A-doubling occurs in spin-forbidden transitions, the inten- 
sity factor for the A -component not shown is obtained by changing the sign of 
all values whose associated transition moment designation contains a 
negative number. 

The three sample cases discussed in this section illustrate the informa- 
tion that is calculated by the computer and printed for both spin-allowed and 
spin-forbidden transitions. The following subsection discusses the range of 
input parameters possible and thereby indicates the comprehensive nature of 
the program. 


Limitations and Capabilities of the Computer Program 

Two types of limitations of the computer program need to be discussed: 

(a) real limitations, for which the computer program does not apply; and 

(b) practical limitations, such as matrix size, that can easily be altered. 

The most important real limitation of the program is the neglect of 
interactions that decouple the orbital angular momentum L from the inter- 
nuclear axis so that the intensity factors for the p- complexes, applicable to 
Hund's case (d) coupling, cannot be determined. Another important real limit 
of the program is that it is only valid for diatomic molecules and for elec- 
tric and magnetic dipole transitions. In fact, the computer program is writ- 
ten for electric dipole transitions but the results also apply to magnetic 
dipole transitions. However, because the parity selection male for magnetic 
dipole radiation is ±*-*±, whereas for electric dipole radiation it is ±*-»-+, 
£ ± -*->E + magnetic dipole transitions must be specified as transitions and 

vice versa. 

There are no important practical limitations to the computer program. 

For example, in nonsigma electronic states, only spin-orbit interactions are 
included in the Hamiltonian, and in sigma electronic states, only the combin- 
ation of spin-spin and second-order spin-orbit interactions are included (see 
appendix C) . However, only if precise wavelengths are desired would a more 
accurate Hamiltonian be necessary. The remaining practical limitations 
actually establish the capabilities of the program. Thus, the computer pro- 
gram is capable of determining intensity factors over the broad range of con- 
ditions listed below. 

1. Maximum spin quantum number of three (i.e., maximum multiplicity of 
seven) . 

2. Maximum of 200 rotation levels in one computer run. 
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3. Spin-allowed transitions with any degree of coupling between Hund's 
case (a) and case (b) . 

4. Spin-forbidden transitions with any degree of coupling between 
Hund's case (c) and case (b) . 

5. Maximum of nine independent transition moments permitted in spin- 
forbidden transitions. 

The maximum multiplicity, the number of rotational levels in a single computer 
run, and the number of independent transition moments permitted are actually 
only limited by the memory size of the computer. Clearly, the capabilities 
of the program are very extensive and permit the calculation of exact inten- 
sity factors for most of the experimentally observed diatomic transitions. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, April 26, 1973 
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APPENDIX A 


LISTING OF COMPUTER PROGRAM 


C THIS PROGRAM COMPUTES INTENSITY FACTORS FOR ALLOWED 

C TRANSITIONS CHONL-LONDON FACTORS ) , AND THE SQUARE ROOT 

C OF THE ROTATIONAL INTENSITY FACTORS FOR FORBIDDEN TRANSITIONS. 

C 

C MAIN PROGRAM 

C 

DIMENSION TITLE1I141 ,TITLE2(121 ,T ITLE3 < 6 1 ,TITLE4( 41 , TITLE5 C 61 , 

1 TITLE6I4) *HEAD2(7,7,21, 

2 HEAD1I71,ALPHAI19,21 
C 

COMMON /CSOLVE/JF.JMIN, J MAX , S I , SF ,L AMI , LAHF, 

1 SIGN1 * ALLOW* EMI SONfSIGSIGt IJ» IHAXI » IMAXF* 

2 NTRANR ( 200 ) » NTRANQ ( 200) * NTRANP (200) » 

3 ITRANI (30*200) » ITRANFf 30*200) • 

4 SAVE<9,i50,200),YH3),YF(3) 

C 

DOUBLE PRECISION SI , SF , L AM I , LAMF, YI , YF , JI , JF, JMI N* JMAX, SI GN1 .SUM 
C 

LOGICAL ALLOW *SIGSIG*HALF«EMISON 
C 

REAL J 
C 

DATA A850R8/1HA/ 

C 

OATA TITLEI/4HS I NG,4HLET ,4HD0UB»4HLET ,4HTRIP,4HLET * 4HQUAR * 

1 4HTET ,4HQUIN,4HTET .4HSEXT.4HET ,4HSEPT,4HET / 

C 

OATA TITLE2/4H SI.4HGMA *4H PI.4H ,4H DE*4HLTA * 

l 4H PH*4HI »4H GA* 4HHMA *4H OT * 4HHER / 

C 

DATA TITLE3/4HEMIS.4HSI0N,4H , 4HABS0* AHRPT I * 4H0N / 

C 

DATA TITLE4/4HINIT.4HIAL *4HFINA*4HL / 

C 

OATA T I TLE5/4HREGU »4HLAR ,4H I NVE , 4HRTED *4H ,4H / 

C 

DATA TITLE6/4HE . U.4HPPER ,4HE. L,4H0WER/ 

C 

DATA HEAD2 

1 /2H1 ,2H12,2Hl3,2Hi4»2Hl5.2H16,2H17 f 

2 2H2 L *2H2 ,2H23, 2H24, 2H25 , 2H26, 2H27, 

3 2H3I *2H32 »2H3 , 2M34* 2H35 , 2H36 .2H37 , 

4 2H4 1 * 2H42 * 2H43* 2H4 , 2H45 , 2H46, 2H47 , 

5 2H51*2H52,2H53,2H54»2H5 ,2H56,2H57, 

6 2H6I »2H62 .2H63* 2H64* 2H65 * 2H6 ,2H67, 

7 2H71,2H72,2H73,2H74,2H7S*2H76,2H7 , 

8 2HI *2H2l »2H3l»2H4l»2H5l *2H61 *2H7I * 

9 2H12.2H2 »2H32*2H42* 2H52*2H62*2H72» 

1 2HL3 » 2H2 3 * 2H3 ,2H43, 2H53, 2H63,2H73, 

2 2HI4»2H24,2H34»2H4 , 2H54.2H64.2H74 . 

3 2H1 5 » 2H25 * 2H35 , 2H45* 2H5 .2H65.2H75* 

4 2H 16 * 2H26 * 2H36 * 2H46 » 2H56* 2H6 ,2H76, 

5 2Hi7*2H27 * 2H37* 2H47 * 2H57 * 2H67 « 2H7 / 

C 

OATA ALPHA/IHZ * IHY * IHX*1HM* IHV , 1HU * 1HT * IHS * 1HR, 1HQ* 

1 IMP * 1H0 * lHNflHM* 1 HL * IHK * 1HJ » LHI * IHH, 1 HH, 

2 LHI • IHJ * IHK * IHL * 1HM, 1HN*1H0*1HP*1HQ*1HR* 

3 IHS* 1HT * 1 HU, IHV * 1HW « 1HX »1HY,IHZ/ 

C 

DATA BLANK/LH /, R/1HR/, Q/1HQ/ * P/1HP / 

C 

C READ IN NEW CASE. 

C 

10 READ! 5*501) SI *LAMI , SIGNI, YI ( II , YII2) ,YI (3) 

REAOI 5*501) SF.LAMF.SIGNF, YF( 1 1 , YF 1 2 > , YF I 3) 

READ! 5*502) JMIN.JMAX 
READ! 5*503) TYPE 
C 

C INITIALIZE LARGE SAVE ARRAY. 

C 

00 15 1=1*200 
DO 15 L=l * l 50 
DO 15 K«l*9 
15 SAVE(K»L»II=0.0 

C 

C INITIALIZE TRANSITION MOMENT ARRAYS. 

C 

DO 17 1=1,200 
NTRANR 111=0 
NTRANQ ( 1 1=0 
NTRANP (I 1=0 
DO 17 L=1 *30 
I TRANI (L * 1 1 = 0 
17 I TRANF1 L *11=0 

C 

C IS THIS AN EMISSION CALCULATION? 

C 

EM I SON 3 . TRUE. 

IF (TYPE.EQ. ABSORB 1 EMISON=. FALSE. 

C 

C IF JMAX IS NOT INPUT, SET JMAX TO JHINUOO. 

C IF JMAX-JMIN IS GREATER THAN 199 SET JMAX TO JMIN*199. 

C 

If UMAX. EQ. 0.000) JMAX=JMIN+100.D0 
DEL JMAX-JMIN 

I F ( DELJ.GT.199. 1 1 JMAX»JMIN*199.0D0 
C 
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C IS THIS A SIGHA TO SIGMA TRANSITION? 

C 

$IG$IG=.TRUE. 

IFI LAM l .GT • .1 .OR .LAKE .GT . • 1 ) SIGSIG 3 . FALSE . 

C 

C IS THIS AN ALLOWED TRANSITION? 

C 

IFUSI-$F).NE..ODO.OR.OABSILAMI-LAHF>.GT.1.1DO) GO TO 20 
C 

IFISIGSIG.AND.SIGNI.NE.SIGNFI GO TO 20 
C 

C THIS IS AN ALLOWED TRANSITION. 

C 

ALLOW 3 . TRUE • 

IMAXI*2.*S!«i.l 
IMAXF*IMAX1 
GO TO 30 
C 

C THIS IS A FORBIDDEN TRANSITION. 

C 

20 ALLOW*.FALSE. 

1HAX1-2.*SI*1.1 

IMAXF»2.*SF+1.I 

C 

C SET SIGN OF SYMMERTRIC MATRIX ELEMENTS IN THE 0 BRANCH MATRIX. 

C 

SIGN1-M.000 

I F C .NOT . S IGSIG) GO TO 30 

IOELS=DABSISF-$!)+.l 

1*1 

IFISIGNI.NE.S IGNF ) 1=0 
SIGNl-<-l)**( IOELSMI 
C 

C IS SPIN HALF INTEGER OR WHOLE INTEGER? 

C 

30 IS-SI+.l 

S1*IS 

HALF*-TRUE. 

I F( DABS! SI-SI ).LT.. 100 » HALF*. FALSE. 

C 

C JMIN ANO JMAX MUST BE HALF INTEGER OR WHOLE INTEGER IN 
C ACCORD WITH THE SPIN. 

C 

JTEST=JM1N*.1D0 

J*JTEST 

IFIHALFI GO TO 35 

I F C DABS I JMIN- J I .GT • • 1 ) JMIN*JHIN«-0.500 
GO TO 37 

35 IF( DABS! JMIN- J I .LT.. 1) JMIN*JMI N*0. 5D0 

37 JTEST»JMAX*.100 

J-JTEST 

IFIHALFI GO TO 39 

I F ( OABSI JMAX-J I .GT..1) JMAX* JMAX+O. 500 
GO TO 40 

39 IF I DABS t JMAX— J I .LT. . 1 1 J MAX=JMAX*0. 5D0 
C 

C FIND MINIHUM VALUES FOR JI AND JF ANO 

C ADJUST JMIN IF NECESSARY. 

C 

40 JI=LAMI— S I 

IF ( JI.GE. O.ODOI GO TO 45 
JI*O.ODO 

IFIHALFI J 1*0.500 
45 IFUJI-JMINI.GT..1) JHIN*JI 

C 

JF*LAMF-SF 

IFIJF.GE. O.ODOI GO TO 50 
JF*0.000 

IFIHALFI JF«0.5D0 
C 

50 IFI JMIN-JFI 55»65,60 

C 

55 IFI IJMIN-JFJ.LT .-1.1001 JMIN* JF-l.ODO 

GO TO 65 
C 

60 IFI IJMIN-JF) .GT. 1.100) JF*JHIN-1 .ODO 

C 

65 CALL SOLVE 

C 

C PRINT HEADING 

C 

WRITE 16*600) 

C 

IFI ALLOW! GO TO 70 
WRITE I 6 t 60l I 
GO TO 75 

70 WRITE 16 (602) 

75 ILAMI*LAMI+.L 

IFIILAMI.GT.5) ILAMI*5 

ILAMF=LAMF+.l 

IFI ILAMF.GT.5I ILAMF=5 

INDEX=4 

IFIEMISON) INDEX*1 
I2=INDEX*2 
C 

INDEX2=5 

IFILAMI.LT. • 100. ANO. SI. LT.. 6 DO I GO TO 77 
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IFISI.LT. .1001 GO TO 77 
I NDEX2* 1 

IFIYI Cl). LT. 0.0001 INDEX2=3 
77 IN0EX3=5 

IFCLAHF.LT. .100. AND. SF.LT..6D01 GO TO 80 

IFCSF.LT. .100) GO TO 80 

INDEX3=1 

IFCYFCH.LT. 0.000) INDEX3*3 
C 

80 WRITEC6.603) T 1TLE5C I NDEX2 ) , T I TLE 5 ( IN0EX2*l)»TITLEl( 2*IMAXI-1 1 

1 TITLE II 2‘IHAXI) ,TITLE2 I 2*1 L AM I *1 ) , T ITLE2 C 2*IL AMI*2 > .SIGN I, 

2 TITLES! INDEX3 I, TITLES ( IN0EX3* 1 ) , T I TLE U 2* l HAXF-1 ) , 

3 TITLE1C2*IHAXF) , T I TLE2 < 2* I L AMF*1 ) ,T ITLE2 C 2* IL AHF*2 ) ,SIGNF, 

4 (TITLE3CI) , I=INOEX, 12) 

C 

IF(IHAXI.EO.l) GO TO 85 
IFC I LAHI.EQ.OI GO TO 82 
WRITE ( 6 « 604 ) TITLE4C11,TITLE4(2),YIC1) 

GO TO 85 
C 

82 IFISI.LT. .900) GO TO 85 
IS=SI 
11-1 

I F I HALF) 11 = 2 

WRITE 16,611) TITLE4C11,TITLE4C2),YII l) 

C 

IFCSI.LT. 1.900) GO TO 85 
12=1 

IF ( SI *GT .2.600) I2»2 
C 

00 84 1=1,12 
13=1*1 

84 WR I TE C6, 6 12 ) I3 ( YI1I*1) 

C 

85 IFC IHAXF.EO.l ) GO TO 90 
IFC ILAMF.EQ.O) GO TO 87 

WRITE ( 6 , 604 ) T I TLE4 C 3 ) ♦ T I TL E4C4 ) , YF C 1 ) 

GO TO 90 
C 

87 IFC SF.LT..900) GO TO 90 
IS-SF 

II = 1 

1 F I HALF) 11=2 

WRITE 16(611) TITLE4C3),TITLE4!4),YFC1) 

C 

IFISF.LT. 1.900) GO TO 90 
12=1 

I F SF .GT .2.600) 12 = 2 
C 

DO 89 1=1,12 
1 3= I* l 

89 WR I TE C6, 612 ) I3,YFCI*1) 

C 

90 IFCLAHI.LT..1 00 .AND. S I .GT . .600) GO TO 95 
IFCLAMF.GT..1D0.0R.SF.LT..6D0) GO TO 100 

95 WRITEC6, 613) 

C 

100 IFCLAHI.LT. .100. OR.LAHF.LT. .100) GO TO 120 
I F ( ALLOW I GO TO 110 
C 

WRITE ( 6,614) 

GO TO 120 
C 

110 WRI TE C6, 615) 

C 

120 INDEX-3 

I F( EM ISON) INDEX = 1 

WRITE (6, 605 1 TI TLE6C IN0EX),TITLE6( INDEX* l ) 

C 

C PRINT DATA 

C 

I I = 2 

I F ( EH ISON ) 11 = 1 
C 

C IS THIS AN ALLOWED TRANSITION? 

C 

I F ( ALLOW) GO TO 190 
C 

C THIS IS A FORBIDDEN TRANSITION. 

C 

IDELS=0 

DELS«SI-SF 

IFCABSCOELSI.LT. .1) GO TO 124 
I F ( OELS.GT .0. t GO TO 122 
IDELS=-l 

IFCOELS.LT. -l.l) I0ELS=-2 
GO TO 124 
122 IDELS=1 

I F COELS.GT. 1. II I D6LS = 2 

C 

C PRINT R BRANCHES. 

C 

124 00 140 K=l, IMAXI 

00 125 1=1, IMAXF 
N 1=9* I DEL S~C K-l )*C I-l)*2M II-l) 

HEA01 CI)=ALPHACNl ,11) 



125 


C 


C 


12B 

130 


C 

140 

C 

C 

C 


145 


C 


c 


148 

150 


C 

160 

C 

c 

c 


165 


C 


C 


168 

170 

C 

180 

C 

C 

c 

c 

c 

c 

c 

c 

190 


195 


C 


IFIHEADK I1.E0.R) HE ADI ( 1) = BLANK 
CONTINUE 

WRITE <6, 6061 I HE ADI I I ) ,R,HEAD2( 1 ,K, II I , I-l, IMAXF » 

DO 140 1 = 1. IJ 
DEL J=I “1 
J=JMIN*D£LJ 
WR l TE ( 6 , 607 ) 

K0=(K-l)*7 

K1-K0*1 

K2=K0*IMAXF 

K3=NTRANR I I I 

DO 130 NTRAN»1,K3 

SUN=O.ODO 

DO 128 KK-KI.K2 

SUM- SUM* SAVE (NTRAN.KK, 11**2 

WRITE (6.608) J, SUM, I TRANFI NTRAN , I ) , ITRANM NT RAN . I ) . 

1 (SAVEINTRAN.KK.n ,KX»Kl,K2) 

CONTINUE 

PRINT 0 BRANCHES. 

DO 160 K-l, IMAXI 
DO 145 I-l, IMAXF 
N1»10*I0ELS-(K-1)*(I-1) 

HEAOl(I)»ALPHA(M,II) 

IFIHEADK n.EQ.O) HE ADI ( I ) = 0LANK 
CONTINUE 

WRI TE (6 . 606 ) ( HE AD 1 (II , Q , HE AD2( I, K , II ) , I-l , I MAXF ) 

DO 160 1-1, IJ 
DELJ-I-1 
J«JMIN*DELJ 
WRI TE ( 6.607 ) 

K0-50*(K-1)*7 

Kl-KO+l 

K2«K0*IMAXF 

K3=NTRANQ ( I I 

DO 150 NTRAN-1.K3 

SUM-0. ODO 

DO 148 KK-K1.K2 

SUM= SUM*SA VE ( NTRAN , KK, I >**2 

WRITE (6,608) J.SUM. ITRANF (NTRAN + 10, 1 ) , I TRAN I ( NTRAN* 10. 1 ) 
1 (SAVEINTRAN.KK.I ),KK=Kl,K2) 


CONTINUE 

PRINT P BRANCHES. 

00 180 K-lt IMAXI 
00 165 I-l, IMAXF 
Nl»ll*IDELS-IK-l)*(I-l)-2*( I I-l | 

HE ADI (I l-ALPHAINl , I I ) 

IFIHEADll D.EQ.P) HEAD1 ( I ) = BLANK 
CONTINUE 

WRITE (6,606) ( HE A01 ( I ) , P , HE AD2( I , K, II) , I-l , IMAXF ) 

00 180 !-l,IJ 
DELJ-I-l 
J-JMIN+DELJ 
WRITE (6,607 ) 

K0-100*(K-1)*7 
K1»K0*1 
K2=K0*I MAXF 
K3-NTRANP ( I ) 

DO 170 NTRAN-1 »K3 

SUM-0.000 

DO 168 KK-KI.K2 

SUM»SUM*S AVE ( NTRAN, KK , I ) **2 

WRITE (6,608) J , SUM, I TR ANF ( NTRAN*20, 1 ) , I TRAN I (NTRAN*20, 1 ) 
1 ( SAVE! NTRAN ,KK , I ) »KK»Kl »K2) 


CONTINUE 

END OF CASE, READ IN NEXT CASE. 

GO TO 10 

THIS IS AN ALLOWED TRANSITION. 

PRINT R BRANCHES. 

DO 210 K-l, IMAXI 
DO 195 I-l, IMAXF 
N 1 =9- ( K- 1 )*( I-l) *2*1 I I-l) 

HEAD1I I)-ALPHA(N1 ,11 ) 

IFIHEADK D.EQ.R) HE ADI ( I ) = BL ANK 
CONTINUE 

WRITE (6,609) (HEADK I ) ,R ,HEA02( I,K,II) , I-l, IMAXF) 

00 210 I-l.IJ 
DEL J= I — 1 


24 



J-JMIN+OELJ 
SUM-0. ODO 
c 

K0»(K-1)*7 

Kl-KO+1 

K2-K0*IHAXI 

DO 200 KK-K1,K2 

K3=KK*50 

K4-KK*100 

ZOO SUM»SUM*SAVE<1»KK,I )*SAVE(l,K3, I ) *S AVE ( 1 , K4 , 1 ) 

WRITE (6,610) J , SUH, ( S AVE ( 1 ,KK , l ) , KK-Kt , K2) 

IF! I ( 1/51*5-1 I.EO.O) WRITE (6,607) 

C 

210 CONTINUE 
C 

C PRINT Q BRANCHES. 

C 

DO 230 K-l, IMAXI 
00 215 l-l, IMAXF 
Nl*10— (K-l !♦( 1-1 ) 

HE ADI ( I ) "ALPHA ( N1 , 1 1 ) 

IFIHEADl (I) .EQ.Q) HE ADI ( I ) -BLANK 
215 CONTINUE 

WRITE (6,609) (HEA01( I), 0, HE AD 2(1, K, II) , I-l, IMAXF) 

C 

00 230 I-l,IJ 
DEL J 3 I- 1 
J-JMIN+OELJ 
SUM-0. 000 
C 

K0=50*( K-l ) *7 

Kl-KO+1 

K2-K0+IMAXI 

DO 22 0 KK-K1 , K2 

K3-KK-50 

K4-KK+50 

220 SUM- SUM* SAVE! 1 ,KK , I )+SAVE(l«K3, I)+SAVE(l»K4,I ) 

WRITE (6, 610) J.SUH, ( S AVEl 1 ,KK , 1 ) , KK-Ki ,K2> 

1F((( 1/5) *5-1 ).EQ.O) WRITE (6,607) 

C 

230 CONTINUE 
C 

C PRINT P BRANCHES 

C 

DO 250 K-l, IMAXI 
OC 235 I-l, IMAXF 
Nl-ll— (K— l)+( I— 1)— 2*( I I-l ) 

HEADK I)-ALPHA(N1,I1 ) 

IF(HEADKI).EQ.P) HEADK I I-BLANK 
235 CONTINUE 

WRITE (6,609 ) ( HE ADI ( I ) , P , HEA02( I,K, II 1 , I-l, IMAXF ) 

C 

DO 250 1=1, IJ 
DEL J 3 I-l 
J-JMIN+OELJ 
SUM-0. ODO 
C 

K 0=1 00* ( K-l I *7 
Kl-KO+1 
K2-K0* I MAX I 
DO 240 KK=K1 ,K2 
K3=KK-100 
K4=KK-50 

240 SUM- SUM* SAVE ( 1 ,KK, I ) ♦SAVE ( 1 »K3 , I ) 4SAVE 1 1 *K4, 1 I 

WRITE (6 ,6 10) J,SUM, (SAVE 1 1 ,KK, I ) , KK-K1 , K2) 

I F ( ( ( 1/5) *5-1 ).EQ. 0 ) WRITE (6,607) 

C 

250 CONTINUE 
C 

C END OF CASE, READ IN NEXT CASE. 

C 

GO TO 10 
C 

C READ FORMAT STATEMENTS. 

C 

501 FORMAT (F3.l,2X,F2.0,2X,Al,3E10.0) 

502 FORMAT (2 F5. 1 ) 

503 FORMAT ( Al ) 

C 

C WRITE FORMAT STATEMENTS. 

C 

600 FORMAT (1H1, 52 X,26H INTENSITY FACTOR PROGRAMME,/ 

1 45X, 41HNAS A-AMES/CR ESS/E E WHITING AUGUST 1972,///) 

601 FORMAT ( 32X,43HSQUARE ROOT OF ROTATIONAL INTENSITY FACTORS, 

1 18H FOR THE FORBIDDEN,//) 

602 FORMAT (48X»35HH0NL~L0ND0N FACTORS FOR THE ALLOWED,/) 

603 FORMAT (27X,2A4, IX,4A4,A1,4H TO ,2A6 , 1 X , * A4 , Al , l 5H TRANSITION IN ,3A4//) 

604 FORMAT ( 35X, 2A4, 47HST ATE SPIN-ORBIT COUPLING CONSTANT - A/B - V = ,F8.3) 

605 FORMAT (/ / ,43X,29HJ VALUE IS FOR THE INITIAL (I,2A4,8HI STATE.) 

606 FORMAT! // ,16X,7HPARTIAL»3X»10HTRANSITION»/IOX,1HJ,7X, 3HSUM ,6X, 6KM0MENT ,7(7X,AI,A1,A2,4X) ) 

607 FORMAT (10X) 

608 FORMAT (7X.0PF 5. I.F11.6.5X, 12 , IH , , 1 1 , 1X,7 1 3X , l PE 12. 5 ) ) 

609 FORMAT (// 10X, 1HJ, 7X,3HSUM,2X, 7(7X, A l,Al ,A2,4XI//) 

610 FORMAT ( 7X,0PF5. 1,F11.6,7(3X,1PE12.5)> 

611 FORMAT (35X,2A4,47HSTAT£ ENERGY SEPARATION « DELTA E!l)/B - ,F8.3) 

612 FORMAT < 75X, 8HDELTA E(,U,6H)/8 » ,F8.3) 



613 FORMAT ( //27X, 58HIN SIGMA STATES WHEN J<S, THE F DESIGNATION FOR THE LEVELS, 

1 /27X.53HNITH THE SAME VALUE OF J ARE ASSIGNEO AS 2S+1 FOR THE, 

2 /27X ,5 1HHIGHEST ENERGY LEVEL, 2S FOR THE NEXT HIGHEST, ETC.) 

C 

614 FORMAT I//27X»59HLANBDA DOUBLING OCCURS IN THE SPECTRUM, BUT THE SQUARE ROOT, 

1 17H OF THE I NT E NS ITY/27X,40HFACTQRS FOR ONLY ONE LINE OF EACH LAMBDA, 

2 40H DOUBLET IS PRINTED. THE VALUES FOR THE/27X , 1 2 HOT HER LAMBDA, 

3 58H COMPONENTS ARE FOUND BY CHANGING THE SIGN OF THE NUMBERS , 

4 IIHPRINTEO, IF/27X,42HTHE TRANSITION MOMENT DESIGNATION CONTAINS, 

5 14H A MINUS SIGN.) 

615 FORMAT I//27X»59HLAMBDA DOUBLING OCCURS IN THE SPECTRUM, BUT THE HONL-LONDON, 

1 21H FACTORS FOR ONLY ONE/27X ,39HL I NE OF EACH LAMBDA DOUBLET IS PRINTED. 

2 39H THE HONL-LONOON FACTORS FOR THE OTHER/27X, 17HLAMBDA COMPONENTS, 

3 14H ARE THE SAME.) 

ENO 
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SUBROUTINE SOLVE 
C 

C THIS PROGRAM: 

C 1-DEVELOPES THE ROTATIONAL HAMILTONIAN. 

C 2-CALLS EIGEN TO GET THE EIGEN VECTORS OF THE HAMILTONIAN. 

C 3-0EVEL0PES THE ROTATIONAL MATRIX. 

C ^-TRANSFORMS THE ROTATIONAL MATRIX. 

C 5-STORES, IN THE LARGE SAVE ARRAY, HONL-LONDON FACTORS 

C FOR ALLOWED TRANSITIONS AND THE SQUARE ROOT OF 

C ROTATIONAL INTENSITY FACTORS FOR FORBIDDEN TRANSITIONS. 

C 

DIMENSION Til*, 14) ,TI I 14, 14) ,TFlt 14,14) , TF2< 14, 14) ,A( 14,14) , 

1 Ctl4»14)»Hf7»7) » SYM (7,7) 

C 

COMMON /CSOLVE/ JF , JMIN, JMAX , S I ,SF,LAH1,LAMF, 

L SIGN I ,ALLOW, EMISON,SIGSIG,IJ,I MAX I , IHAXF , 

2 NTRANRI200 ) , NTRANQ ( 200 ) , NTRANP ( 200) , 

3 IT RANI (30,200) ,ITRANFI30,200), 

4 SAVEI9, 150, 200), YI<3),YF(3> 

C 

CGMHON/CE I GEN/HH ( 28) ,TT(49) ,N 

C 

OOUBLE PRECISION A, DEL J , JMI N, JMAX , J , J I , JF , Y( 3 ) , Y I , YF , LAM, LAMI , 

1 LAMF,S,SI,SF,H,HH,W,HOLD,TT,T,TI ,TFl ,TF2 , OMEGA, OHEGAI , C , 

2 OMEGAF , E I , EK , SUM, FACTOR , SIGN, SIGN!, SIGN 2, SIZ£,SYM, 

3 Cl 

LOGICAL ALLOW , EMI SON, SIGSIG, FIR ST 
C 

REAL INITAL 
C 

DATA R/1HR/ , Q/1H0/ • P/1HP/, F1NAL/1HF/ , INITAL/lHI/ 

C 

INTEGER SAVE1 , SAVE 2 , SAVE3 , SAVE4 , SAVE 5, SAVE 6 
C 

C INITIALIZE MATRICES. 

C AND SET FLAG FOR FIRST COMPUTATION. 

C 

DO 10 1*1,14 

00 10 K* l » 14 
T I ( K i I ) =0.000 
TF1IK, 11*0.000 
TF2(K,ri=0.000 
C IK , l ) *0.000 

10 A ( K , I ) =0 . ODO 
FIRST*. TRUE. 

S AVE3*0 
NSAVE-0 
C 

C START PRIMARY CALCULATION AT THIS POINT. 

C THE INTENSITY FACTORS ARE CALCULATED IN 

C TERMS OF THE INITIAL VALUE OF J, IE. IN TERMS OF 

C J-LOWER FOR ABSORPTION AND IN TERMS OF 

C J-UPPER FOR EMISSION. 

C 

DO 900 I J = 1 ,200 
DEL J= I J- 1 
J I- JMI N+OEL J 
C 

C SET VALUES FOR INITIAL LEVEL CALCULATION. 

C 

J = JI 
LAM*LAMI 
S = SI 

1 MAX* IMAX I 
STATE* INITAL 
DO 15 1*1,3 

15 YII)»Yim 

GO TO 30 
C 

C SET VALUES FOR FINAL LEVEL CALCULATION. 

C 

20 J=JF 

LAM=L AMF 
S=SF 

I MAX* I MAXF 
STATE=F1NAL 
DO 25 1*1,3 
25 YIII-YFm 

c 

C INITIALIZE HAMILTONIAN AND TRANSFORMATION MATRICES. 

C 

30 DO 35 1*1,7 

DO 35 K=l ,7 
35 HIK, 11*0.000 

DO 40 1=1,14 
DO 40 K=1 , 14 
40 TIK ,1 )*O.ODO 

00 45 1*1,28 
45 HHII >=0.000 

DO 47 1=1,49 
47 TTII)»O.ODO 

C 

C THE HAMILTONIAN MATRIX IS PEDUCEO AT LOW J VALUES. 

C SET INDECIES DEFINING WHICH ROWS AND COLUMNS HAVE NDN-ZERO ELEMENT 

C 

OMEGA=LAM+S 

IFI IJ+.ll.LT. OMEGA I GO TO 50 
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) 


c 

C HAMILTONIAN IS FULLY DEVELOP! D. 

C 

11 = 1 
12*1 MAX 
GO TO 110 
C 

C HAMILTONIAN IS REDUCED. 

C 

50 T£ST=OMEGA 

00 60 1*1 *IMAX 
1 1*1*1 

TEST-TEST-1.0 

IF! ( J ♦ - 1 > •GT.TEST) GO TO 70 
60 CONTINUE 

11= I MAX 
GO TO 80 
C 

70 IFU.GT..1) GO TO 100 

C 

C THERE IS ONLY ONE ELEMENT IN THE HAMILTONIAN. 

C SET TRANSFORMATION MATRIX. 

C 

80 12=11 

C 

90 Tll,ll» <-l)**Il 

N*1 

GO TO 2*0 

C 

C THERE may BE MORE THAN ONE ELEMENT IN THE 
C HAMILTONIAN MATRIX. 

C 

100 I2»2.*J*.l 

1 2= 11+12 

IFd2.GT.IMAX) 12 = 1 MAX 

C 

c IS THERE MORE THAN ONE ELEMENT IN THE 

C HAMILTONIAN MATRIX? 

C 

110 IFU2.EQ.il) GO TO 90 
C 

C THERE IS MORE THAN ONE ELEMENT IN THE 

C HAMILTONIAN MATRIX. SET SIZE OF MATRIX. 

C 

N«I2-Il*l 

C 

C DEVELOP UPPER HALF OF SYMMETRICAL HAMILTONIAN MATRIX. 

C J(J-H) AND SISU) TERMS IN THE 01 AGONAL ELEMENTS ARE NOT INCLUOED. 

C 

K=0 

13=12-1 

DO 120 1=11,13 
K»K*1 
E I * I 

H<K.K«1)— OSORTI < J-OMEGA+EI)*(J*OMEGA-em .000 )*F I * I 2. ODO*S-E I ♦! . 000 ) ) 
HI K+l , K)=H(K »K + 1 1 

120 HIK ,K)=-<OMEGA-Eld.DO)**2-CS-Em.DO)**2 
C 

E 1*12 
K=K+1 

HlK.K) =-l OMEGA-E 1*1 .001**2- I S-E1*1.D01 **2 
C 

C IS THIS A SIGMA STATE? 

C 

K»0 

IFILAM.LT. .100) GO TO 140 

C 

C THIS IS NOT A SIGMA STATE. 

C INCLUOE SPIN-ORBIT INTERACTION ( Y=A/B) . 

C 

DO 130 1=11,12 
K=K+1 
E I® I 

130 H(K,K)=H(K,K)+Y( l I *LAM*( S-EI+1.00) 

GO TO 152 
C 

C THIS IS A SIGMA STATE. 

C INCLUDE SPIN-ORBIT ♦ SPIN-SPIN ENERGY SEPERATI ONS IDELTA E/B). 

C 

140 DO 150 1=11,12 
K=K* 1 

C 

GO T0I152, 152, 141, 142, 143,144, 145), IMAX 
C 

141 IFU.EQ.21 GO TO 150 
13=1 

GO TO 146 
C 

142 1FII .E0.2.0R. I .EQ.3) GO TO 150 
13-1 

GO TO 146 
C 

143 IFII.E0.3) GO TO 150 
13=2 

IF( I.EQ.2.0R.I.E0.4) 13=1 
GO TO 146 
C 



144 IFI I.EQ.3.0R.I.C0.4I GO TO 150 
13*2 

IF ( 1 .E0.2.0R. I .E0.5) 13=1 
GO TO 146 

145 IFI1.EQ.4I GO TO 150 
13*3 

IFI I.NE.l.OR.I .NE.7) 13=2 
IFI I.EQ.3.0R.I.EQ.51 13*1 

C 

146 HIK,K1=HIK,K1*Y< 131 
150 CONTINUE 

C 

C IS SYMMETRIZING MATRIX NECESSARY? 

C NEEOED FOR NON-SINGLET SIGMA STATES IN FORBIODEN TRANSITIONS. 

C 

152 I F ( ALLOW. OR. LAM. GT ..1.0R.S.LT..1) GO TO 225 

C 

C CONSTRUCT SYMMETRIZING MATRIX. 

C 

IFIN.EO.NSAVE) GO TO 180 
NSAVE-N 
00 155 1-1,7 
00 155 K-1,7 
155 SYM( I »K ) =0. 000 
C 

Cl *.7071 0678 100 
K=0 
C 

160 K*K*1 

KK-N+l-K 

I F ( K.EQ.KK) GO TO 170 
SYMI K ,K1 »Cl 
SYMIKK , K ) *C 1 
SYMI K « KK) =C1 
SYMIKK, KK)*-Cl 
C 

IFUKK-Kl.GT.il GO TO 160 
GO TO 180 
C 

170 S YM(K , K1 =1 .000 
C 

C SYMMETRIZE HAMILTONIAN MATRIX IF REQUIRED. 

C 

180 DO 200 1=1, N 
DO 200 K-l»N 
C ( K , I ) “0 . ODO 
00 190 KK*1,N 

190 CIK»I1«C(K,II*HIK,KK)*SYHIKK,I) 

200 ’ CONTINUE 
C 

DO 220 1=1, N 
DO 220 K=1,N 
HI K , I 1 = 0.000 
00 210 KK=l,N 

210 HIK,I)-HtK,l )*SYNIK,KK)*CIKK, I 1 
220 CONTINUE 
C 

C TRANSFORM UPPER HALF OF SYMMETRICAL HAMILTONIAN INTO A ONE -DIMENSIONAL 

C ARRAY, HH TO BE COMPATIBLE WITH EIGENFUNCTION SUBROUTINE. 

C 

225 KK=0 

DO 230 I-1,N 

00 230 K= 1,1 
KK-KKU 

230 HH( KK I =H( K, 1 1 
C 

C CALI SUBROUTINE TO FIND THE El GENFUNCTI DNS OF 
C THE HAMILTONIAN MATRIX BY JACOBI-VON NEUMAN METHOD. 

C EIGENFUNCTIONS ARE ORDERED IN A ONE-DIMENSIONAL ARRAY, TT WITH 

C THE LARGEST EIGENVALUE FIRST, ETC. 

C 

CALL EIGEN 
C 

C EXPANO THE ON E-OIMEN SIGNAL ARRAY, TT INTO THE TWO-DIMENSIONAL 

C TRANSFORMATION MATRIX, T THAT 01 AGONAL I Z ES THE HAMILTONIAN. 

C ENSURE CONSISTANT PHASE OF THE EIGENFUNCTIONS BY MAKING 

C THE FIRST ELEMENT IN EACH COLUMN NEGATIVE. 

C 

KK»0 

DO 235 1=1, N 
C 1*1 .000 

IFITTIKKM1.GT. 0.) Cl — 1. ODD 

DO 235 K»t,N 

KK-KK+1 

235 T(K,I 1=C1*TT(KK1 
C 

C IF THIS IS A FORBIDDEN TRANSITION AND IF 

C THIS STATE IS NOT A SIGMA STATE, EXPAND THE 

C TRANSFORMATION MATRIX TO ALLOW FOR LAMBDA DOUBLING. 

C 

240 NN*N 

I F ( AL LOW .OR . LAH.LT • • 1 1 GO TO 260 

C 

C EXPAND THE TRANSFORMATION MATRIX. 

C 

NN=2*N 

1 I-NN+1 
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00 250 1=1, N 
DO 250 K*1 , N 
K3=II-K 
13*11-1 

250 TIK3,I3)*T(K, I ) 

C 

C TRANSFER T MATRIX INTO INITIAL ( T I > OR FINAL (TF) 

C MATRIX. SET FLAGS AND INDECIES FOR PROPER BRANCHES 

C PRELIMINARY TO DEVELOPING THE ROTATIONAL MATRIX. 

C 

C IS THIS THE INITIAL OR THE FINAL STATE? 

C 

260 IFtSTATE.EO. FINAL) GO TO 280 
C 

C THIS IS THE INITIAL STATE. 

C TRANSFER THE T MATRIX INTO THE T! MATRIX. 

C SET INDECIES FOR INITIAL STATE. 

C 

00 270 1*1, NN 
00 270 K*1 , NN 
270 Til t , Kl *T ( l yKI 
C 

111*11 

112*12 

N1*N 

C 

C IF THIS IS THE FIRST TIME THRU OR IF THE 0 BRANCHES 

C HAVE NOT BEEN COMPUTED, GO DIRECTLY TO THE FINAL STATE. 

C 

IF (FIRST ) GO TO 20 
IF(BRANCH.Nc.Q) GO TO 300 
C 

GO TO 330 
C 

C THIS IS THE FINAL STATE. 

C ARE CONDITIONS SET FOR 0 BRANCHES? 

C 

280 IFIOABSI JI-JFKGT..1 ) GO TO 310 
C 

C CONDITIONS ARE SET FOR 0 BRANCHES. 

C INVERT T MATRIX AND PUT INTO TF1 MATRIX. 

C 

00 290 1*1, NN 
DO 290 K*1 , NN 
290 TFl (1 ,K) *T(K, I ) 

C 

C SAVE 1ND1CIES DENOTING SIZE OF HAMILTONIAN AND 

C WHICH ELEMENTS ARE FILLED. 

C 

SAVEi*Il 
SAVE2*I2 
SAV£3*N 
FIRST*. FALSE. 

C 

C INCREASE FINAL STATE ROTATIONAL QUANTUM NUMBER 

C BY 1 ANO COMPUTE THE NEW FINAL STATE. 

C 

300 JF*JF+l.ODO 
GO TO 20 
C 

C THIS IS EITHER A P OR AN R BRANCH. 

C NORMALLY JF WILL BE GREATER THAN Jl AT THIS POINT. 

C HOWEVER, THE FIRST TIME THRU Jl MAY BE GREATER. 

C 

310 IF(JF.GT.JI) GO TO 350 
C 

C FIRST TIME THRU. INVERT T MATRIX AND PUT INTO TF1 

C MATRIX AND SET INDICIES. 

C 

Do 320 1=1, NN 
DO 320 K= 1 , NN 
320 TF1 ( I , K) =T( K« I ) 

C 

IF1*II 
I F2* 1 2 
NF=N 
C 

C IF THIS IS AN EMISSION CALCULATION, WE ARE READY 

C TO COMPUTE R BRANCHES. IF ASSORT ION, P BRANCHES. 

C 

330 I F( EM ISON) GO TO 370 
C 

C THIS IS A P RRANCH CALCULATION. 

C 

3A0 BR ANCH*P 
C 

C SET INDEX FOR SAVING P BRANCHES IN THE LARGE- SAVE ARRAY 

C 

I K*100 

S IGN2=- 1 • 000 
GO TO 390 
C 

C INVERT T MATRIX, PUT INTO TF 2 MATRIX ANO SET INDICIES. 

C 

350 DQ 360 1*1, NN 
OC 360 K=1,NN 
360 TF2(1,K)=T(K,I) 
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c 

IP1*I 1 
IF2-I2 
NF*N 
C 

C IF THIS IS AN EMISSION CALCULATION* WE ARE READY 

C TO COMPUTE P BRANCHES. IF ABSORPTION, R BRANCHES. 

C 

I F I EM ISON I GO TO 340 

C 

C THIS IS AN R BRANCH CALCULATION. 

C 

370 BRANCH=R 
C 

C SET INDEX FOR SAVING 0 BRANCHES IN THE LARGE SAVE ARRAY. 

C 

IK*0 

SIGN2 — 1.0D0 
GO TO 390 
C 

C THIS IS A 0 BRANCH CALCULATION. 

C 

380 BRANCH=0 
I K=50 

SIGN2=*l.0D0 

C 

C DEVELOP ROTATIONAL MATRIX. 

C 

C SE.T INDICIES AND SIGN OF LAMBDA. 

C 

390 1=0 

E 1 = 1 1 2 
NI I = N I 
NFF=NF 
NTRAN-0 
SIGN® 1 . ODO 
C 

395 I F f ALLOW. OR. LAMP. LT..1) GO TO 4.00 

S I GN=- l • 000 
NFF=2*NF 
C 

C SET INOICIES AND OMEGA FOR THE INITIAL STATE, 

C COLUMNS OF THE MATRIX. 

C 

400 1=1*1 

E I = EI — l.ODO 
OMEGA I=LAMI+SI-E I 
1 1 =N I -I + 1 
C 

K =0 
EK=IF2 
C 

C SET INDICIES AND OMEGA FOP. THE FINAL STATE, 

C ROWS OF THE MATRIX. 

C 

410 K=K+I 

EK=EK- l . OCO 
KK=NFF-K+1 

I F < SIGN. LT.O. > EK= (KK-K-ll/2 
OMEGAF=SIGN*LAMF*SF-EK 
C 

C FIND OELTA OMEGA ANO DELTA LAMBDA. 

C 

DELO=OMEGAF-OMEGA I 
DELL-LAHF-LAMI 
C 

C FOR ALLOUEO TRANSITIONS DELTA OMEGA MUST EQUAL DELTA LAMBDA. 

C 

IF (.NOT. ALLOW I GO TO 415 
I F ( AB5( DELO-DE LL ) .GT •• l ) GO TO 510 
C 

C WHAT TYPE OF BRANCH? 

C 

415 IF ( BRANCH. £0.0) GO TO 470 

IF( (JF-JI ).GT. .11 GO TO 440 
C 

C P BRANCH IN ABSORPTION OR 

C R BRANCH IN EMISSION. 

C 

I F ( J I .LT..I) GO TO 810 
I F ( ABS(OELO) -LT..11G0 TO 430 
I F ( OELO.GT . 0. ) GO TO 420 
I F ( 0EL0.LT.-1.1JG0 TO 510 
C 

Af KK * 1 1 1 =-DSQRT ( ( J I IOMEGA I ) *( J1 +OMEGAI - 1 • 000 1 / 1 2. 000* J II) 

GO TO 500 
C 

420 IF(OELO.GT.l.l) GO TO 520 

C 

A I KK, 1 1 ) =DSQRT ( ( J I-OMEGA I I *( J I -OMEGA I -1 . ODOI / ( 2. ODO* J I ) 1 
GO TO 500 
C 

430 A (KK, 1 1 ) =DSORT( ( J1 ♦OMEGA I ) * ( J I -OMEGA I I / J 1 1 
GO TO 500 
C 

C R BRANCH IN ABSORPTION OR 

C p BRANCH IN EMISSION. 
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c 

440 I F C ABS(DELO) .LT..11G0 TO 460 
TFCOELQ.GT.O.I GO TO 450 
Z F < DELO.LT.-l.l )G0 TO 510 
C 

A(KK,I! >=DSQRT(< J I-OHEGA I ♦ 1 .000 1*1 J! -OMEGA 1*2. ODO) / 1 2. ODO* ( J 1*1 .0001 )1 
GO TO 500 
C 

450 IF( OELO. GT .1.1) GO TO 520 
C 

A ( KK, ! I I c-oSQRTC ( J I ♦OMEGA I + 1 .ODO ) * ( J I ♦OMEGA I *2 .0001/ < 2.OD0* ( JI ♦ 1 . 000 ) ) ) 
GO TO 500 
C 

460 AIKK.II )=>DSQRT( I J! ♦OMEGA 1*1 .000 ) *( J I-OHEGA I ♦ l.ODO) /< J I ♦ l . 000) ) 

GO TO 500 
C 

C 0 BRANCH IN ABSORPTION OR EMISSION. 

C 

470 IFUI.LT..1) GO TO BIO 

1 F ( ABS ( OELO) . LT . .1 ) GO TO 490 
l F ( DELO.GT. 0. ) GO TO 480 
IF ( OELO. LT. -I. I) GO TO 510 
C 

A ( KKt I ! ) a DSQRT ( ( J I +OMEGA I ) * ( J I -OMEGA 1*1. 000) * (2.0D0*JT+1. 000) 

1 /( 2.000* JI* (JI+1.000))) 

GO TO 500 
C 

480 I F ( DELO.GT. 1.1) GO TO 520 

C 

A(KK»II)=DSQRT( (JI~OMEGAI)*(J I ♦OMEGA I+l.0D0)*(2.0D0*JI+l. ODO) 

1 /<2.000*Jl*( JI+1.000))) 

GO TO 500 
C 

490 A IKK, II ) “0HEGA1 *DSQRT((2. 000* J 1*1 .ODO )/(JI*(JI+l .000) ) ) 

C 

C IF THIS IS AN ALLOWED TRANSITION LOOP BACK UNTIL ALL 

C ELEMENTS ARE IN THE ROTATIONAL MATRIX. IF THIS IS 
C A FQRBIODEN TRANSITION, FIND THE SYMMETRICAL ROTATIONAL 

C MATRIX ELEMENT AND TRANSFORM THEM TOGETHER. 

C 

500 IF( .NOT. ALLOW) GO TO 540 

C 

510 IF(K.LT.NF) GO TO 410 
IFIK.GT.NF) GO TO 525 
C 

520 1 F ( SIGN .GT • .0) GO To 530 

C 

C IS FORBIDDEN BRANCH COMPLETEO? 

C 

I F ILAMI .GT. .1.0R.!I+1).LT.I1) GO TO 522 
IFCI.EQ.II) GO TO 810 
C 

522 S IGN»*l .ODO 
EKMF2 
GO TO 410 
C 

525 IF(K.LT.NFF) GO TO 410 
C 

530 IF(I.EQ.NI) GO TO 535 
IF ( ALLOW) GO TO 395 
C 

C IS FORBIDDEN BRANCH COMPLETED? 

C 

IF(LAMI.GT..1.0R.( I + D.LT.II) GO TO 395 
GO TO. 810 
C 

C IF THIS IS AN ALLOWED TRANSITION, MUST NOW TRANSFORM 

C THE ROTATIONAL MATRIX. IF THIS IS A FORBIDDEN 

C TRANSITION, THIS BRANCH IS COMPLETE. 

C 

535 I F l ALLOW ) GO TO 580 

GO TO 810 
C 

C THIS IS A FORBIDDEN TRANSITION. FIND THE SYMETRICAL 

C ELEMENT WITH PROPER PHASE FACTORS. CHECK SIGMA TO SIGMA 

C TRANSITION FOR EXISTENCE OF TRANSITION MOMENT WHEN BOTH 

C OMEGAS EOUAL ZERO. 

C 

540 IF(SIGSIG) GO TO 560 
C 

IFILAMI.LT. .1) GO TO 565 
C 

C INITIAL STATE IS NOT A SIGMA STATE. 

C 

N1!«2*NI 
I 3=N! ♦ ! 

GO TO 567 
C 

C THIS IS A SIGMA TO SIGMA TRANSITION. 

C ARE BOTH OMEGAS EQUAL TO ZERO? 

C 

560 IF! DABS (OMEGA I) .LT..l.AND.DABSf OMEGAF ).LT..l) GO TO 570 
C 

C AT LEAST ONE OMEGA IS UNEQUAL TO ZERO, SET SYMMETRICAL ELEMENT. 

C 

565 13=1 

567 A(K,I3)«SIGN1*SIGN2*A(KK, I I) 



GO TO 571 


C 

C BOTH OMEGAS ARE ZERO. 

C THERE IS NO SYMMETRICAL ELEMENT. 

C FOR THE TRANSITION MOMENT TO EXIST SIGNl MUST BE NFGATIVE. 

C 

570 I F (SIGNL . LT .0 . ) GO TO 580 

C 

C THE TRANSITION MOMENT CORRESPONDING TO THIS MATRIX 

C ELEMENT DOES NOT EXIST. He CAN ALLOW FOR THIS BY 

C SETTING THE MATRIX ELEMENT EQUAL TO ZERO. ALSO THIS 

C BRANCH IS COMPLETE. 

C 

AIKK, 1 1 1*0. ODO 
GO TO 810 
C 

C TRANSFORM THE ROTATIONAL MATRIX. 

C 

C SYMMETRIZE FORBIDDEN TRANSITIONS. 

C 

571 IFIK.EG.KK) GO TO 576 
IF( I3.EQ.III GO TO 575 
1FIK.GT.KK) GO TO 573 
IF( I3.GT.il) GO TO 572 

C 

HH<1)»A(K,I3)*0.5D0 
HHl 2) *0.000 
HH( 3 )*0.000 
HH(4)*AIKK, 111*0.500 
GO TO 578 
C 

572 HHI1)*0.0D0 

HHl 21 °-A IK, 131*0.500 
HHI3)»A(KK,II)*0.5D0 
HH14)*0.0D0 
GO TO 578 
C 

573 I F I I3.GT.II ) GO TO 574 
HH ( l ) =0.000 
HH(2)*A|KK, 111*0.500 
HH(3)=-AIK, 131*0.500 
HHI4) =0.0D0 

GO TO 578 
C 

574 HHl l ) *A( KK, II 1*0.500 
HHI 2) -O.ODO 

HH 1 3 1 *0. ODO 
HHI4)=AIK,I3)*0.5D0 
GO TO 578 
C 

575 HH(1)*A(K,I3>*0. 70710678100 
HHI2I = AI KK, 13 1*0. 70710670100 
HHl 3) =0. 000 

HH 1 4 ) ”0.000 
GO TO 578 
C 

576 I FI I3.GT.tl) GO TO 577 
HHU) -AIK, I 31*0.707106781 DO 
HHl 2) =A|K,1 I) *0.70710678100 
HHl 3) >0.000 

HH(4) »0. 000 
GO TO 578 
C 

577 HH(l) »A(K, II )*0. 70710678100 
HHI2)=-A IK, 13 >*0.70710678100 
HHl 3 ) “0.000 

HHl 4) *0.000 
C 

578 A(K,I3)*HHI D*HHI2)*HHI3) *HH l 4 ) 

AlKK*Il) * HHl 1 ) -HHl 2 ) — HHl 3>*HH|4) 

IFIK.EQ.KK.0R.13.EQ.il) GO TO 580 
A IK, I I )*HH(1)-HHI2)*HHI3)“HHI4) 
A(KK,I3)»HHI1)*HHI2)-HH(3)-HHI4) 

C 

C MULTIPLY ON THE RIGHT BY Tt. 

C 

580 DO 600 1 1*1 • NX I 
DO 600 K 1 *1 ,NFF 
CIK1,I1)*O.ODO 
DO 590 12-1, Nil 

590 CIKI,II)*C(K1,I1)*AIK1, I2)*TI ( 12,11) 

600 CONTINUE 
C 

C MULTIPLY ON THE LEFT BY TF 

C 

IFIBRANCH.EQ.Q.OR.I JF-JI I.LT..I) GO TO 630 
C 

C FOR P BRANCHES IN EMISSION AND R BRANCHES IN 

C ABSORPTION MULTIPLY ON LEFT BY TF2. 

C 

DO 620 Il-l.NII 
DO 620 Kl-l.NFF 
A(K1,I1)*0. 000 
Do 610 12*1, NFF 

610 A(Kl,m-A<KI,lll+TF2(Kl, 121 *Ct 12,11) 

620 CONTINUE 
GO TO 655 
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c 

c 

c 

c 

630 


640 

650 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

655 


C 


660 

c 

662 


663 
C 

664 


665 
C 

C 

C 

C 

C 

C 

666 

C 

C 

C 


C 

C 

C 


670 

C 

680 

C 

650 

C 

695 

C 

C 

C 


FOR 0 BRANCHES AND FOR P BRANCHES IN ABSORPTION 
AND R BRANCHES IN EMISSION MULTIPLY ON LEFT BY TFl. 

DO 650 11=1, Nil 
00 650 Kl=l,NFF 
A (Kl, 1 1) =0.000 
DO 640 12=1, NFF 

A(Kl«lL}=A(Kl»!L)+TFL(Klf 12) *C ( 12,11) 

CONTINUE 

FOR ALLOWED TRANSITIONS LAMBDA DOUBLING WAS NOT INCLUDED. THEREFORE , THE 
MATRIX ELEMENTS FOR SIGHA-PI TRANSITIONS MUST BE MULTIPLIED BY THE SQUARE 
ROOT OF 2. 

FOR FORBIOOEN TRANSITIONS; 

IF ONLY THE INITIAL STATE IS A SIGMA STATE, THE RESULTS 
IN THE UPPER AND LOWER HALVES OF THE MATRIX MUST BE COMBINED. 

IF ONLY THE FINAL STATE IS A SIGMA STATE, THE RESULTS IN 
THE LEFT AND RIGHT HALVES OF THE MATRIX MUST BE COMBINED. 

IF NEITHER STATE IS A SIGMA STATE, THE RESULTS 

ARE EITHER IN THE UPPER-LEFT AND LOWER-RIGHT QUADRANTS, DR IN 

THE UPPER-RIGHT AND LOWER-LEFT QUADRANTS. THEREFORE, WE CAN ALWAYS GET 

ONE LAMBDA COMPONENT BY AODING THE MATRIX ELEMENTS IN THE UPPER-LE^T 

AND LOWEP-LEFT OUADRANTS I ONE OF WHICH IS ZERO). THE OTHER LAMBDA 

COMPONENT CAN 86 FOUND BY SIMPLY CHANGING THE SIGN OF ALL TERMS 

WHOSE TRANSITION MOMENT IS DESIGNATED WITH A MINUS SIGN IN THE 

FIRST LOCATION, IE I-X,Y>. 

IF(SIGSIG) GO TO 666 
IFI. NOT. ALLOW) GO TO 662 

IFUAMI. GT.. IDO. AND. LAMP. GT.. IDO) GO TO 666 

Cl»l. 41421356200 
00 660 K3=i,NF 
00 660 13=1, N! 

A(K3,I3)=Cl*A(K3,I3) 

GO TO 666 

1FUAHF.LT. .100) GO TO 664 
DC 663 13=1, N! 

DO 663 K3=l ,NF 
K4=NFF+I-K3 

AIK3,I3)=AIK3,I3)+A ( K4» 13) 

GO TO 666 

00 665 13=1, NI 
00 665 K3=1,NF 
I4=NI 1*1-13 

A(K3»I3)=A(K3»I3)*A(K3,I4) 

FOR STATES OTHER THAN SIGMA STATES, IF THE HAMILTONIAN 
IS NOT FULLY DEVELOPED (SMALL J VALUES) THEN THE ELEMENTS 

IN THE TRANSFORMED ROTATIONAL MATRIX ARE SHIFTED TO 

ENABLE PROPER LABELING OURING OUTPUT. 

NDELI-0 

IFUAMI. LT..1. OR. III. EQ.l) GO TO 695 
DETERMINE COLUMN SHIFT. 

NOEL I = 111-1 

I F f V 1 C I I .LT.O.ODO) NOEL 1 = I MAX I- 1 Il-N!*l 
IF(NDELI-LE.O) GO TO 695 

SHIFT COLUMNS 

DO 690 Kl=l ,NF 

00 670 11=1, NI 
I2=NI*l-I 1 

1 3*1 2+NDEL I 

A(Kl , 1 3 ) = A ( Kl , I 2 ) 

00 680 11=1, NOEL I 
A (K 1 , 1 1 1=0.000 

CONTINUE 

NDELF=0 

IFUAMF.LT..I .OR. IF1.EQ.1 ) GO TO 735 
OETERHINE ROW SHIFT. 

NDELF-IFl-l 

IF( YFI 1) .LT.O.ODO) NOELF® I MAXF- 1 Fl-NF+l 
I F ( NDELF.LE .0 ) GO TO 735 

SHIFT ROWS. 

DO 730 11=1, IMAX! 

DO 710 Kl=l ,NF 
I2=NF+1-K1 
13=1 2+NDELF 
A(I3,U) = A( 12,11) 

00 720 Kl =1 , NDELF 
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720 A(Kl,HI=O.ODO 
C 

730 CONTINUE 
C 

C STORE THE TRANSFORMED ROTATIONAL MATRIX ELEMENTS 

C AND INITIALIZE THE A MATRIX FOR A NEW CALCULATION. 

C 

735 K4=N0ELF +NF + 1 

I 4=NDEL I+NI+1 
K5= I MAXF-NOELF-NF 
1 5= IMAX I-NOEL I-NI 
C 

C IS THIS AN ALLOWED TRANSITION? 

C 

IFI. NOT. ALLOW! GO TO 760 
C 

C THIS IS AN ALLOWED TRANSITION. 

C STORE THE SQUARE OF THE A MATRIX ELEMENTS. 

C THESE ARE THE HONL-LONDON FACTORS FOR A SINGLE LAMBDA COMPONENT. 

C 

DO 755 I I s 1 1 N I 
16*14-! 1 
DO 750 K1=I »NF 
K7»(I5»!i-l )*7+ 1K*K5+K 1 
K6=K4— K1 

750 SAVE! I , K7, I J I = A { K6 f 1 6 ! **2 

755 CONTINUE 

DO 757 11=1,14 
DO 757 K i=l » 14 
757 AtKl, 111*0.000 
C 

GO TO 810 
C 

C THIS IS A FORBI ODEN TRANSITION. 

C RECORD THE NUMBER OF ELECTRONIC TRANSITION 

C MOMENTS AND THEIR DESIGNATIONS. 

C 

760 NTRAN=NTRAN+1 

ie=ei+.i 

KE=EK*.l 
K6=NTRAN*IK/5 
I TRAN 1 1 K6 , I J ) = 1 MAX I- 1 c 
C 

I F { LAMP. LT ..II GO TO 780 
IF( S IGN.LT.O. ) GO TO 770 
C 

I TRANF IK6,I J1=*K— NF+IMAXF-IF2 
GO TO 790 
C 

770 ITRANFIK6, I Jl *K- IMAXF* I F 1-2 
GO TO 790 

C 

780 ITRANF1K6»IJ>=!HAXF-KE 

C 

C IS THIS AN EMISSION CALCULATION? 

C 

790 I F C .NOT • EMISON) GO TO 799 
C 

C THIS IS AN EMISSION CALCULATION. AOJUST DESIGNATION 

C OF THE TRANSITION MOMENT TO CORRESPOND TO ABSORPTION. 

C 

K 1 =♦ 1 

I F ( SIGSIG I GO TO 794 
!F(LAMF.LT. .11 GO TO 792 
lFILAMI.LT.. 1> GO TO 791 

C 

IFISIGN.LT. 0.) K 1 *— 1 
GO TO 796 
C 

791 IF ( SIGN .GT.O. I GO TO 796 
C 

U-ITRANI IK6, IJ) 

I TRAN M K6 ,IJ»*-I TRANF CK6, 1 J ) 

I TRANF (K6, 1 J) -I MAXI + 1-11 
GO TO 799 
C 

792 IFIOMEGAF.LT. -.1) GO TO 796 
C 

I 1*ITRANF(K6, IJI 
I TRANF (K6 , I J ) *-ITRANI IK & , I J I 
I TRAN I (K6,IJ)*I MAXF+1- 1 1 
GO TO 799 
C 

794 I F ( OMEGA F . GT . . 1 ) GO TO 798 
C 

796 I 1 = ITRANF <K6» 1 J I 

I TRANF (K6, 1 J)=K1*ITRANI (K6.IJI 
I IRANI ( K6 ? I J I *K 1 *1 L 
GO TO 799 
C 

798 1 1 = I TRANI ( K6 , I J ) 

I TRAN I (K6,IJ) = IMAXFM-1TRANF(K6,IJJ 
I TRANF (K6, 1 J I = I MAXI ♦ l- 1 1 

C 

C SAVE A MATRIX ELEMENTS. 

C 

799 DC 000 I i*I »NI 
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CO 800 Kl = l ,NF 

K7=U5*I1-1I*7*IK*K5*K1 

K6=K4-Kl 

16=14-11 

800 SAVEtNTRAN.KT, IJI=A(*6, 16) 

C 

DO B05 ! 1*1 * 14 

00 805 KI*i ,14 

805 A(K1«I1)*0.0D0 

C 

GO TO 510 
C 

C A SET OF BRANCHES HAS BEEN COMPUTEO. 

C STORE NUMBER OF ELECTRONIC TRANSITIONS MOMENTS 

C IN COMPUTATION JUST FINISHED AND SET FLAGS * 

C IDICIES* ETC APPROPRIATE FOR NEXT SET OF BRANCHES. 

C 

810 IFI8RANCH.NE.0) GO TO 820 
C 

C A SET OF 0 BRANCHES HAVE BEEN COMPUTED. 

C NEXT STEP IS TO INCREASE JI AND COMPUTE A NEW INITIAL STATE. 

C 

NTRANQ1I J)*NTRAN 
GO TO 890 
C 

820 lF(BRANCH.EO.P) GO TO 830 
C 

C A SET OF R BRANCHES HAVE BEEN COMPUTED. 

C 

NTRANRU J)*NTRAN 

1 F ( EMISON ) GO TO 840 
GO TO 860 

C 

C A SET OF P 8RANCHES HAVE BEEN COMPUTED. 

C 

830 NTRANP ( I J I = NTRAN 

I F ( EMISON ) GO TO 860 
C 

C REPLACE TFI MATRIX WITH TF2 MATRIX BEFORE INCREASING JF 
C AND COMPUTING NEW FINAL STATE. 

C 

040 K4» SAVE3 

I F I .NOT . ALLOW. AND. LAM. F *GT .. 1 ) K4»2*$AVE3 
CO 850 11=1, K4 
DO 850 K1*1,K4 
850 TFKKltl I )*TF2f Kl ,IU 
C 

JF* JF+1. 000 
GO TO 20 
C 

C IS THIS THE FIRST CALCULATION? 

C 

860 IFCFIRSTI GO TO 870 
C 

C THIS IS NOT THE FIRST CALCULATION. SAVE INDICIES 

C AND COMPUTE 0 BRANCHES NEXT. 

C 

SAVE4* I F 1 
SAVE5-IF2 
SAVE6-NF 
IF1-SAV61 
IF2-SAVE2 
NF=SAVE3 
SAVE1 *S AVE4 
SAVE2-SAVE5 
SAVE3«SAVe6 
GO TO 380 
C 

C THIS IS THE FIRST CALCULATION. SAVE INDICIES ANO 
C TFZ MATRIX ANO COMPUTE NEW INITIAL STATE. 

C 

870 FIRST*. FALSE. 

SAVE1» JF1 
SAVE2*l?2 
SAVE3=NF 
C 

DO 080 11*1, NFF 
DO 800 Kl=i,NFF 
880 TFl(Kl f Il)*TF2(Kl f Il) 

C 

C WOULD NEXT VALUE OF JI 8E GPEATER THAN J MAX? 

C 

890 IF( ( JMAX-JI) .LT..9) GO TO 910 

C 

C END OF BIG 00 LOOP 

C 

900 CONTINUE 
C 

I J*200 
910 RETURN 
END 



SUBROUTINE EIGEN 


PURPOSE 

COMPUTE EIGENVALUES AND EIGENVECTORS OF A REAL SYMMETRIC 
MATRIX 

USAGE 

DESCRIPTION OF PARAMETERS 

A - ORIGINAL MATRIX (SYMMETRIC), DESTROYED IN COMPUTATION. 
RESULTANT EIGENVALUES ARE DEVELOPED IN DIAGONAL OF 
MATRIX A IN DESCENDING ORDER. 

R - RESULTANT MATRIX OF EIGENVECTORS (STORED COLUMNWISE, 

IN SAME SEQUENCE AS EIGENVALUES) 

N - ORDER OF MATRICES A AND R 

REMARKS 

ORIGINAL MATRIX A MUST BE REAL SYMMETRIC (STORAGE MOOE-l) 
MATRIX A CANNOT BE IN THE SAME LOCATION AS MATRIX R 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIREO 
NONE 

METHOD 

DIAGONALI ZATION METHOD ORIGINATED BY JACOBI AND ADAPTED 
BY VON NEUMANN FOR LARGE COMPUTERS AS FOUND IN 'MATHEMATICAL 
METHOOS FOR OIGITAL COMPUTERS', EDITED BY A. RALSTON AND 
H.S. WILF, JOHN WILEY AND SONS, NEW YORK, 1962, CHAPTER 7 


SUBROUTINE EIGEN 
CCMMON/CE IGEN/A(28),R(49) ,N 
C 

DOUBLE PRECISION A,R, ANORM, ANRMX,THR,X,Y, SINX,S!NX2 ,COSX, 
1 COSX2,SINCS, RANGE 

C 

C GENERATE IDENTITY MATRIX 

C 

5 RANGE»1.0D-12 

10 IQ»-N 

00 20 J«l,N 
IQ= IQ+N 
DO 20 1-1, N 
1J-I0*! 

P(I J)=O.ODO 
IF(I-J) 20,15,20 
15 R(IJI~1.0D0 
20 CONTINUE 
C 

C COMPUTE INITIAL AND FINAL NORMS (ANORM AND ANORMX) 

C 

25 ANORH-O.ODO 
00 35 !-l,N 
00 35 J-I.N 
IF(I-J) 30,35,30 
30 I A= !♦ ( J*J- J )/2 

ANCRN=AN0RM»A( IA)*A( IA) 

35 CONTINUE 

I F ( ANORM) 165,165,40 
40 AN0RM=DSQRT(2.000*AN3RMI 

ANRMX=ANORM*RANGE/FLO AT ( N ) 

C 

C INITIALIZE INDICATORS ANO COMPUTE THRESHOLD, THR 

C 

IND-0 

THR-ANORM 

45 THR«THR/FLOAT(N) 

50 L»l 
55 M«L*1 
C 

C COMPUTE SIN ANO COS 

C 

60 M0=(M*M-M)/2 
L0»(L*l-L)/2 
LM-L*NQ 

62 IF(DABS(A(LM))-THR) 130,65,65 
65 INO-1 
LL-L+LQ 
MM-M+MQ 

X=0.5D0*(A(LL)-A(MM)) 

69 Y — A(LM1/DS0RT( A(LM)*A(LM)*X»X) 

1 F ( X I 70,75,75 

70 Y=- Y 

75 SINX«Y/PS0RT(2.000*( 1 .000*- ( DSQRT ( 1.000-Y*Y)) ) ) 

SINX2»SINX«SINX 
78 COSX-OSQRTI 1.OD0-SINX2) 

cosx2»cnsx*cosx 

SINCS -SI NX*COSX 

c 

C ROTATE L AND M COLUMNS 

C 

I LQ-N* ( L — 1 1 
I MQ-N*( M— l ) 

DO 125 1*1, N 
1 Q* ( I* I-I )/2 
IF(I-L) 30,120,80 
80 IF(I-M) 85,120,90 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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65 IMM+NO 
on 70 «J5 
-50 1 M= M*- 1 0 

95 IFIt-LI 100,105.105 
100 I L= I + LO 
GO TO 110 
105 IL*L*!0 

110 X = AUL)*CO$X-A( IMJ*SINX 

A( 1 H) =AC IL)*SINX*AUM>*COSX 
AIIll*X 
120 ILR® ILQM 
I MR = I MO* I 

X=R(ILR)*C0SX-R( IMR)*SINX 
R ( I MR 1 = R C ILR 1 *$ INX*R( IMR»*C0SX 
R(ILR)«X 
125 CONTINUE 

X=2.0D0*A(IM) *SINCS 
Y = ACL L)*C0SX2«-A ( MN)*S INX2-X 
X*A ( LL I *S INX2+A ( MM) *C 0SX2+X 

A I LMI = ( A ( LL I- A( MM) )*SINCS+A(LM)*( C0SX2-S INX2 ) 

AILLI-Y 

A|MM)»X 

C 

C TESTS FOR COMPLETION 

C 

C TEST FOR H « LAST COLUMN 

C 

130 IF(M-N) 135.140.135 
135 M=M*1 

GO TO 60 
C 

C TEST FOR L = SECOND FROM LAST COLUMN 

C 

140 IFIL-IN-lll 145.150.145 
145 L=L*l 

GO TO 55 

150 1 F ( I NO- 1 ) 160,155,160 
155 1 ND=0 

GO TO 50 
C 

C COMPARE THRESHOLD WITH FINAL NORM 

C 

160 I F ( THR- ANRMXI 165,165,45 

C 

C SORT EIGENVALUES AND EIGENVECTORS 

C 

165 I 0=-N 

00 165 1*1, N 
IQ*IQ+N 

LL- !♦ « I*I-I»/2 
J0*N*{ 1-2) 

DO 165 J*I,N 
J0= JO+N 
MM=J* ( J* J-J 1/2 

1 F l A ( LL J - A l MM 1 ) 170,185,185 
170 X=A (LL I 

A ( LL ) =A ( MM) 

A(MM)=X 

175 00 180 K-l.N 
ILfi-IQ+K 
IMR=JQ+K 
X=R( ILRJ 
R ( ILR 1 =R ( IMR I 
180 R ( I MR I =X 
185 CONTINUE 
RETURN 
END 
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APPENDIX B 


DETAILED FLOW CHARTS OF COMPUTER PROGRAM 


PERTINENT BRANCHING SYMBOLS 





> BRANCHING ON SAME PAGE 


/ 



TO PAGE 


> BRANCHING BETWEEN PAGES 



FROM PAGE 


/ 































































APPENDIX C 


INTERACTION TERMS INCLUDED IN HAMILTONIAN MATRIX 


In nonsigma states the spin-orbit interaction is of first order and is 
usually the dominant interaction. In sigma states the spin-orbit interaction 
is very small and, in fact, is on the same order as the spin-spin interactions 
Therefore, nonsigma and sigma electronic states are discussed separately. 


NONSIGMA ELECTRONIC STATES 


The diagonal elements, of the Hamiltonian matrix in terms of the 

basis functions for nonsigma electronic states, are given by 

= BJ{J + 1) - BQ 2 + 55(5 + 1) - BZ 2 + AA.Z. (Cl) 

'Z' 'V % 'V 

where A is the spin-orbit coupling constant. The result is the same for both 
lambda substates; therefore, for this discussion A. can be considered positive 
Equation (Cl) can be written L 

A = W * 1) + S(S * 1) - (A + ZJ 2 - Z\ + lkZ i (C2) 

where Y = A/B, or as 

ff . Cy) = j(j + 1) + 5(5 + 1) - A 2 - 2Z 2 + (J - 2) AE . (C3) 

Bn^ i i 

The criterion for regular or inverted electronic states is specified by 
the sign of A or J; that is, if A > 0, the state is regular, and if A < 0, the 
state is inverted. 


SIGMA ELECTRONIC STATES 


The spin-orbit interaction of sigma electronic states is of second order 
and A = 0. Therefore, the spin-orbit and spin-spin interactions are of the 
same order and both should be considered. Both of these interactions depend 
on the | £ | , and their combined effect is included in the following analysis. 


The diagonal elements of the Hamiltonian matrix in terms of the basis 
functions for sigma electronic states can be written 




5(5 + 1) - 


„ A E. 

2E 2 + -=£■ 
% B 


(C4) 
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where AE^ is the separation of the spin substates for N = 0 caused by spin- 
spin and second-order spin-orbit interactions. Figure 14 illustrates this 
separation and indicates the degeneracy of the rotational levels with the same 
value of N and |e|. The number of energy separations is obviously related to 
the spin multiplicity. There are, for example, zero energy separations for 
singlet and doublet states, one for triplet and quartet states, two for quar- 
tet and sextet states, etc. 

The values of A E- in equation (C4) are found by extrapolating the energy 
separations shown in figure 14, to N = 0. The reference level for AE^ is 
usually chosen to be the matrix element(s) with the lowest value of |z|. 

Whether a sigma electronic state is a regular or an inverted state is 
determined by the sign of A E. If A E > 0, the state is regular, and if AE < 0, 
the state is inverted. 
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TABLE 1.- ROTATIONAL MATRIX ELEMENTS 


[The nonvanishing matrix elements <fiV where <a> = f (J' ;J,Q) . 

The factors for a given matrix element are taken from different rows of the 
same column of this table. The choice of columns depends on the value of J’-J. 
In all cases, the first factor is taken from row one and the second factor from 
row two or three.] 


Factor 

R Branch 

Q Branch 

P Branch 

J* = J + l 

j 1 - j 

j ' = j - 1 

f (j\ j) 

r t i/2 

[4(J + I)j 

r n-i/2 

[4J(J+I)/(2J+I)j 

r t-i/2 
[4J] 

g(J', fi; Jil) 

r l 

2[(J + il + l)( J-il + l)J 

zn 

r t 1/2 

g(j' , ft±l ; jn) 

+ [2(j±n+i)(j±a + 2)] l/2 

[2(j + n)(j±n+i)] l/2 

±[2(j + a)(j + ii-i)] l/2 


TABLE 2.-<R a> BASIS MATRIX FOR A SPIN-FORBIDDEN 2 n «- 4 E + TRANSITION 
e 


\ 4 z + £ 

\ CS 

\ 14 

, \ v> 
2 n \_s. 

<A s' j'l \ 

/X 

roicvj 

roicv) 

rOlM 

+ 

O 

■<~s 

:?<* 

IOICVJ 

O 

-Km 

& 

POICVJ 

+ 

O 

^X 

"S 

micvi 

rOICVI 

micM 

O 


Re, “1 

R e 2 a 2 



®-mA 


Re 4 a 4 

R e 5 a5 


<-'iiK-44 


R e 6 “6 

R e 7 a 7 

R e 8 “8 




R e 9 a 9 

R e| 0 a IO 


TABLE 3.- <R 0 a> P - BRANCH BASIS MATRIX FOR A SPIN-FORBIDDEN 2 II 4 I + TRANSITION 


r 7T 

\ 4 y + “"2 

\ L a 

\. ! W 
o \ y> 
2 n \ji 

/X 

~~3 

rO|C\l 

rolcv) 

rO|OJ 

+ 

Q 

A 

—loo 

rOlOJ 

+ 

O 

A 

“3 

—leu 

4 

i 

POICVJ 

-f 

_o_ 

<x 

"~3 

PO/CVJ 

1 

rO|CO 

1 

fO|c\j 

+ 

O 

turn 

^e| a l 

R e 2 a 2 



<-HK*4 

R e 3 Ct 3 

R e 4 °4 

R e 5 °5 




_R e 5 a 5 


- R e 3 a3 




- R e 2 c>Z 

’Re, “1 
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TABLE 4.- SIMPLIFIED FLOW CHART OF COMPUTER PROGRAM 














TABLE 5.- A DIAGRAMMATIC REPRESENTATION OF THE ELECTRONIC-ROTATIONAL 
HAMILTONIAN FOR A 2 n ELECTRONIC STATE IN TERMS OF BASIS 
FUNCTIONS < AST \<Qj\ 


\ A 

\ 

\ C5 

\ 

X w 

X c/) 

\ — 

<AS2|<fij|\ 

-x 

io|cj 

/X 

—ICO 

—ICO 

—ICO 

—ICO 

—ICO 

* 

—ICO 

rO|fM 

-ICO 

<'HKM 

H|, 

H I2 



<«HKW 

h I2 

h 22 






H 22 

H|2 




H|2 

H,l 


TABLE 6 .-</? a> BASIS MATRIX FOR THE BRANCHES OF A 2 IR->- 4 II TRANSITION 
e 


\ A 

\ 4n s 

X X 

o \ h 

2 n \ 

X— 

<A s' 2 '|<n'j 1 \ 

/X 

“3 

m|f\j 

7 \ 

fOI<\i 

(?) '?- 0J 

<rfl<!i'l 

© 

/X 

“3 

—|CJ 

— ICO 
1 

^2^fOlco 

< r f-Kf-| | l 

© 

X 

— ICJ 

/X 

rOlOO 

rOICO 

i 

“3 

-ICO 

1 

7\ 

-ICO 

rOKvi 

X 

rO|CO 

i 

/X 

— ICO 
rOlOO 

X 

U3|CJ 

/X 

rO|CJ 

rotco 


R e (2,4) 

xa(2,4) 

R e (2, 3) 
xa(2, 3) 

R e (2, 2) 
xa!2, 2) 


R e (-2, i) 
xa(-2,l) 




<'r?Ki J IO 


R e (l, 3) 
xa(l,3) 

R e (l, 2) 
xa(l,2) 

R e (l. 1) 
xa(l,l) 

R e H,l) 

xa(-l,l) 

R e (-I,2) 
x a(-l,2) 






R«(-l,2) 

xal-1,3 

R e (-1, 1) 
xa(-l,l) 

R e (l, 1) 
xa(l,l) 

R e 0, 2) 
xa(l,2) 

R e ( 1 , 3) 
xa(l,3) 






R e <-2, 1) 
xa(-2,l) 


R e (2, 2) 
xa(2,2) 

R e (2, 3) 
xa(2,3) 

R e (2,4) 

xa(2,4) 





TABLE 7.-<R a> BASIS MATRIX FOR THE P-BRANCHES OF A 4 Z + -*-^ 4 Z" TRANSITION 
e 


4 Z“ C 

\ b 

4 ^ 

\ y> 

<A'S , Z , kfl'J'| 

'S- 

rO>CM 

rOK'J 

rOjOJ 

®o_ 

—IC\J 

/\ 
— IOJ 
rO|CJ 

© Jo. 

— <C\J 
1 

/X 
Ho d 

fO|OJ 

©_o_ 

POICM 

• 

7\ 

rO|OJ 

1 

rO|CVJ 

©k 

<o + |||<|j-i|© 

-R e (!.D 

x a(l,l) 

-R e (l, 2) 
x a(l, 2) 




-R e <2,l) 
xa( 2,1) 

-R e !2,2) 
x a(2,2) 

R e (3,2) 
x a(3,2) 


<o*| 4|<4 j - | |© 


-R e (3,2) 

xa(3,2) 

R e (2, 2) 
x a(2,2) 

R e (2, 1) 

x a (2,1) 

(o’i-iK-iH© 



Red, 2) 
x a(l, 2) 

ReO.D 
x a(l,l) 


TABLE 8.- 


DESCRIPTION OF INPUT CARDS FOR PROGRAM 


Card 

number 




Number of 


Columns 

Format 

Content 

cards per 

Comments 




case 



1-3 

F3.I 

Spin 

2 

The first card contains 


6-7 

F20 

A 


the information for the 


10 

A! 

+ or - ; i.e. symbol for 


initial state and the 

■ 



Z states 


second card contains 


11-20 

EIO.O 

Y for nonsigma states. 


the information for the 

■ 



AE(I)/B for Z states. 


final state. 

1 

21-30 

EI0.0 

AE(2)/B forZ stated 




31-40 

EIO.O 

AE(3)/8 for Z states. 






*as needed 





S|gfl 


1 








3 


D 

A for absorption. 

1 

It is recommended 




Anything foremission. 


that ABSORPTION or 



■ 



EMISSION always be 
used. 
















TABLE 9.- COMPUTER PRINTOUT FOR SAMPLE CASE 1 


HHNL-LONDC1N FACTORS FOR THE ALLOWED 

INVERTED TRIPLET SIGMA ♦. tq REGULAR TRIPLET PI TRANSITION IN EMISSION 


INITIAL STATE ENERGY SEPARATION = DELTA E(!>/5 = -10.000 

FINAL STATE SPIN-ORBIT COUPLING CONSTANT = A/B = Y * 100.000 


IN SIGMA STATES WHEN J<$, THE F DESIGNATION FOR THE LEVELS 
WITH THE SAME VALUE OF J ARE ASSIGNED AS 2S*1 FOR THE 
HIGHEST ENERGY LEVELt 2S FOR THE NEXT HIGHEST f ETC. 



J VALUE .IS 

FOR THE INITIAL 

(IE. UPPER) 

STATE. 

J 

SUM 

R! 

0R12 

PP.13 

0.0 

0.000000 

0.00000 

0.00000 

0.00000 

1.0 

*5 . 999999 

9.52267E-01 

0.00000 

0.00000 

?.o 

9.99959? 

1.345946 00 

9. 770885-02 

0.00000 

3.0 

13.99999? 

1.70604E 00 

3.11 3B6E-01 

2.58351E-01 

4.0 

1 7.999998 

2.05673E 00 

6.09323P-01 

5.28580E-01 


J 

SUM 

SR21 

P.2 

OP 23 

0.0 

0.000000 

0.00000 

0.00000 

0.00000 

1.0 

5.999999 

l.OOOOOE 00 

0.00000 

0.00000 

?.0 

9.999998 

1.4 994 OE 00 

5 • 99281 F-04 

0.00000 

3.0 

13.999997 

1.9S679E 00 

4. 34552E-0? 

3.72I96E-01 

*.0 

17.999998 

2.4 9098 E 00 

1.381 19F-02 

7.45205E-01 

J 

SUM 

TR3 1 

SR32 

P 3 

0.0 

2.000000 

0.00000 

0.00000 

0.00000 

1 .0 

5.999999 

4 . 77330F-02 

0.00000 

o.oocoo 

2.0 

9.999999 

1 • 53B60E-01 

9 . 0249 l E-OT 

0.00000 

3.0 

13.999999 

2.94775E-01 

1.685515 00 

7.72769E-02 

4.0 

17.999998 

4.47505C-0I 

2.37*40^ 00 

2.294ft ! E-OI 


J 

SUM 

01 


P012 

no? 2 


0.0 

0.000000 

0.00000 


0.00000 

0.00000 


1.0 

5.999999 

1.44596E 

00 

1.25636P-01 

0.00000 


2.0 

9.9999Q7 

2.29541 p 

00 

5 • 390225-01 

1.4261 5 e 

00 

3.0 

13.599997 

3.08029E 

00 

1 .1 79 3 7c 00 

2 . 2 5773E 

00 

4.0 

1 7.999998 

3.84303E 

00 

1.956975 00 

2 

00 


J 

SUM 

R02^ 


02 

P022 


0.0 

0.000000 

0.00000 


0.00000 

0.00000 


1.0 

5.999999 

1 .49940E 

00 

5.99281F-04 

0.00000 


2.0 

9.999999 

2.49P97E 

00 

8 • 59462F-0? 

1.663105 

00 

3.0 

13.999997 

3.48453E 

00 

2 . 9 1 786E-02 

2.90295E 

00 

4.0 

t 7.999998 

4.46575E 

00 

6.70398F-02 

4.017?1 E 

00 


J 

SUM 

S031 

P.022 


03 

0.0 

2.000000 

0.00000 

0.00000 - 


0.00000 

1.0 

5.999999 

5.4fc360E-02 

2 . 87 3765 

00 

0.00000 

2.0 

9.999999 

2.0561 8E-01 

4.45094P 

00 

2.45525F.-01 

3.0 

l ? . 9999°9 

4.35160E-01 

5.7809 35 

00 

6.7 k ! 73E-01 

4.0 

17.999998 

6.912155-01 

6.97250F 

00 

1.23663F 00 


J 

SUM 

PI 

0P12 

NP 1 3 


0.0 

0.000000 

0.00000 

0.00000 

0.00000 


1.0 

5 . 999999 

4.94976F-0! 

1.63153 c -0l 

2.R1800E 

00 

?.o 

9.999997 

9.51723E-01 

4 . 953695-01 

2 • 64 867E 

00 

3.0 

13.999997 

1.37756E 00 

9.313115-01 

2.89796E 

00 

4.0 

17.999998 

1.79069E 00 

1.415B3F 00 

2.94918E 

00 


J 

SUM 

0P21 

P2 

OP23 


0.0 

0.000000 

0.00000 

0.00000 

0.00000 


1.0 

5.999999 

4.981 81F-01 

5 .693985-03 

2.9961 2E 

00 

2.0 

9.999998 

9.93554E-01 

1 .789875-02 

3.32! 89E 

00 

3.0 

13. 099997 

1.43495E 00 

3 .909415-02 

? • 72596E 

00 

4.0 

17.909998 

1.97123E 00 

7 . 1 6605F-02 

4.1571 IE 

00 


J 

SUM 

RP31 

QP32 


P? 

0.0 

2.000000 

7.99041E-04 

1.99920E 

00 

0.00000 

1.0 

5.999999 

9.2380RE-03 

2.83136E 

00 

1 . 832 68E-01 

2.0 

9.999999 

S.95046E-02 

3.48773E 

00 

4.90331 E-01 

3.0 

13.999999 

1.45440E-01 

4.03135E 

00 

8.66373E-01 

4.0 

1 7.599998 

2 .499706-01 

4.51507E 

00 

1.27925E 00 
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TABLE 10.- PARTIAL PRINTOUT FOR SAMPLE CASE 2 


SQUARE ROOT OF ROTATIONAL INTENSITY c *CT0RS FOR THE FORBIDDEN 
REGULAR QUINTET SIGMA - TO INVERTED TRIPLE 1 PI TRANSITION TV tPSDRPTION 


INITIAL STATE ENFPGY SEPARATION = DELTA Ell)/? 

DELTA El?)/? 

FINAL STATE SPIN-ORBIT COUPLING CONSTANT = A/P = Y 


s.ocn 

? O.OPT 
-SO. 000 


IN SIGMA STATES WHEN J<S, THE F DESIGNATION FOP. THE LEVELS 
WITH THE SAME VALUE OF J ARE ASSIGNED AS 2S*i FOP THE 
HIGHEST ENERGY LEVEL, ?S FOR THE NEXT HIGHEST, ETC. 



J 

VALUE IS 

FDR THE INITIAL 

(IE. LOWER) STA 

T -« 


PARTIAL 

TRANSITION 



J 

SUM 

MOM F NT 

SP1 

TR2J. 

UR 3 1 

0.0 

0.000000 

-2 ,3 

0.00000 

0.00000 

0.00000 

0.0 

0.000000 

-1,3 

0.00000 

0.00000 

0.00000 

1.0 

0.000000 

-3,2 

0. 00000 

0.00000 

0.00000 

1.0 

0.000000 

-2,2 

0. 00000 

0.00000 

0.00000 

1.0 

0.000000 

-1,2 

0. 00000 

0.00000 

0.00000 

1.0 

0.000000 

1.2 

0. 00000 

0.00000 

0.00000 

t.o 

0.000000 

-2.3 

0. 00000 

0.00000 

0. 00000 

I .0 

0.000000 

-1 ,3 

D. 00000 

0.0000° 

0.00000 

2.0 

0.117959 

-3, I 

-2.422R3E-01 

-2.826925-0? 

-1.4’549f-03 

2.0 

0.023592 

-2.1 

1.26445E-02 

-1 . 52 34 C E-0I 

-1 . 40931 F-02 

2.0 

1.423921 

-3.2 

-1.18922F DO 

-9. 921791-0? 

-4 . 96 ?£ Q F -0 3 

2.0 

1-139137 

-2.2 

-8.78636E-02 

1.05363F 00 

1.03*69! -01 

2.0 

0.4271 77 

-1.2. 

2.525565-03 

-6.3381 4E-02 

6 . 505 02 F-01 

2.0 

0.4271 77 

1.2 

-2 . 5 2556E-03 

6.33P14F-02 

-6.505025-01 
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Figure 3.— Unitary transformations that transform the Hamiltonian matrix from basis functions 
to symmetrized basis functions for 3 2 and 4 2 electronic states. 
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Figure 4.— Designation of the rotational levels for a 4 n electronic state. 
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Figure 5.— Designation of the rotational levels for a 3 2 electronic state. 
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Figure 6.— Letter designations for changes in N and J during a transition. 
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Figure 7.- Illustration of branch designation. 
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Figure 8.— Branch designation scheme for spin-allowed 4 I1 •«-> 4 2 transitions. 
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Figure 9.— Branch designation scheme for spin-forbidden 2 fl «-»■ 4 fl transitions. 
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Figure 10.— Possible nonzero matrix elements in the transformed rotational matrix of a 3 2 -*■ 3 X 

transition; X any nonsigma electronic state. 
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Figure 1 1.— Possible nonzero matrix elements in the transformed matrix of a 3 X-> 3 Y transition; 

X and Y any nonsigma electronic states. 


F 7 F 6 F 5 F 4 F 3 F 2 F l 


F 7 

E 3 

1 

Iriverte 

d 

1 

1 

1 


F 6 


E 2 j 


l 

l 

1 


F 5 


: E ' 


1 

Regular 


f 4 


i 

i 


1 

1 

1 


F 3 


~ i 
i 
i 


1 

1 

l 


F 2 


i 

i 

i 


1 

l 


F l 


i 

i 

i 


1 

l 

1 



(a) 


+ A Hamiltonian submatrix for the J = l rotational levels of a 


7 I1 electronic state 


F l F 2 F 3 F 4 F 5 
= - 2 -10 12 
N = 7 J = 10 9 8 7 6 

6 1 8 7 6 5 

5 8 7 6 5 4 

4 Z 6 5 4 3 

3 § 5 4 3 2 

2 5 4 3 2 

I 4 3 2 UE 3 

3, £ UEg 

2 | I^Ej 0 


F 6 

3 
5 

4 
3 


(b) Designation of the rotational levels 
for a regular electronic state. 
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Figure 1 2.— Shifting of the J = 1 matrix elements of a 7 n electronic state so that the position of 
the elements in the matrix correspond to the designations for the fully developed matrix. 
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Figure 13.— Sketch of the SAVE (9, 150, 200) array. 
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Figure 14.— Illustration of spin splitting in sigma electronic states. 
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